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Abstract 

New  techniques  for  calculating  goodness-of-fit  statistics  for  normal  distribu¬ 
tions  with  parameters  estimated  from  the  sample  are  investigated.  Samples  are 
generated  for  a  Normal(0,l)  distribution.  Critical  vadues  are  calculated  for  five  mod¬ 
ifications  to  the  Anderson-Darling  statistic  and  five  modifications  to  the  Cramer- Von 
Mises  statistic.  An  extensive  power  study  is  done  to  test  the  power  of  the  new  statis¬ 
tics  versus  the  power  of  the  unmodified  statistics. 

Powers  of  six  of  the  new  statistics  show  minimal  to  no  improvement,  two  of 
the  new  statistics  show  a  marked  decrease  in  power,  aind  two  of  the  new  statistics 
show  an  overall  increase  in  power  over  the  unmodified  statistics.  One  of  these  two 
improved  statistics  was  the  obvious  better  of  the  two,  and  it  was  a  modification  to 
the  Anderson- Darling  statistic. 

Complete  tables  of  critical  values  for  sample  sizes  n=4  through  n=53  are  in¬ 
cluded  for  all  Anderson- Darling  and  Cramer- Von  Mises  statistics,  both  modified  and 
unmodified. 
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MODIFIED  ANDERSON-DARLING  AND  CRAMER- VON  MSES 
GOODNESS-OF-FIT  TESTS  FOR  THE  NORMAL 
DISTRIBUTION 

1.  Introduction 

1.1  bachgrov.T..d 

The  amount  of  money  the  Air  Force  receives  from  Congress  to  accomplish  the 
mission  continues  to  siirink.  With  technological  advances,  problem  solving  often 
becomes  more  and  more  complex.  It  is  not  enough  to  just  solve  the  problems,  it  is 
essential  that  they  be  solved  as  efficiently  and  effectively  as  possible.  When  these 
problems  involve  quantifiable  data  on  which  analysis  can  be  accomplished,  decision¬ 
makers  rely  on  analysts  to  assist  them.  One  problem  solving  technique  often  used 
by  an«Jysts  to  help  decision-makers  is  simulation.  Simulation  is  a  modelling  of  the 
real-world  that  incorporates  the  critical  characteristics  of  the  real-world  problem. 
Simulation  is  much  less  expensive  than  conducting  live  tests  and  can  be  modified  as 
the  problem  changes.  Decision-makers  use  the  results  of  the  analyst’s  simulation  to 
make  efficient  and  effective  decisions. 

Simulation  models  are  made  up  of  many  different  features  and  parts.  One  of 
these  features  is  that  they  mimic  at  least  one  real-world  random  phenomena.  For 
example,  the  random  arrival  of  customers  at  an  Air  Force  Accounting  and  Finance 
office  would  be  a  feature  in  a  simulation  used  to  examine  customer  service  times. 
In  order  for  the  simulation  to  properly  mimic  the  real-world  random  phenomena, 
a  sample  of  real-world  data  must  be  collected  and  its  characteristics  scrutinized. 
(21:43).  After  collecting  the  data,  the  analysts  compare  this  data  to  a  theorized  cu¬ 
mulative  distribution  function,  i.e.  a  sample  is  estimated  to  be  from  the  exponential 


1-1 


distribution  with  -  10.  At  this  point,  the  analysts  want  to  know  how  well  the 
data  fits  the  hypothesized  distribution.  To  check  if  the  data  fits  the  hypothesized 
distribution,  a  test  called  a  goodness-of-fit  test  (GOFT)  is  performed.  GOFTs  are  a 
class  of  statistical  comparisons  between  the  chciracteristics  of  a  sample  of  data  and 
the  characteristics  of  a  estimated  theoretical  distribution  with  the  hope  that  the  two 
match. 

A  GOFT  is  used  to  find  the  appropriate  distribution  of  the  sample  data  so 
that  assumptions  and  predictions  can  be  made  about  the  behavior  of  the  population. 
This,  in  turn  allows  claims  to  be  made  about  the  population  and  subsequently,  about 
the  validity  of  the  simulation.  For  example,  "With  95  percent  confidence,  I  predict 
this  missile  will  destroy  the  target.” 

A  GOFT  is  conducted  by  going  through  these  four  steps. 

1.  Collect  a  Scimple  of  data  for  the  event  of  interest. 

2.  Estimate  a  distribution  of  this  sample  data. 

3.  Choose  the  type  of  GOFT  to  use  and  compute  the  appropriate  test  statistic. 

4.  Compare  this  test  statistic  to  the  critical  value  from  statistical  tables. 

if  the  test  statistic  is  less  than  or  equal  to  the  table  value,  then  the  estimated 
distribution  is  not  rejected.  If  the  test  statistic  is  greater  than  the  table  value, 
then  the  estimated  distribution  is  rejected  and  the  GOFT  is  repeated  with  another 
estimate  of  the  distribution.  When  a  test  shows  a  good  fit,  it  gives  a  degree  of 
confidence  to  the  analyst  that  the  hypothesized  distribution  describes  the  real-world 
data  and  that  it  can  be  used  effectively'  in  a  simulation  model  or  any  other  analysis. 
Accuracy  is  critical  at  this  stage  of  the  game;  no  matter  how  sophisticated  the 
simulation  model  is,  the  choice  of  an  incorrect  distribution  for  a  random  process  can 
make  the  ensuing  analysis  useless  (17) 

There  is  not  just  one  test  that  will  tell  the  analyst  whether  or  not  an  estimated 
distribution  is  correct.  Over  the  yeais,  many  different  ty'pes  of  GOFTs  have  been 
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developed.  With  the  development  of  each  new  GOF'T,  the  goal  was  to  create  a 
GOFT  that  Wtis  better  than  the  hkst.  The  most  commonly  used  GOFTs  are  the 
Chi-square,  the  Kolmogorov-Smirnov  (KS),  the  Anderson-Darling  (A-D)  and  the 
Cramer-von  Mines  (C-VM)  (7:109-114).  Many  of  the  more  recent  developments  with 
GOFTs  have  been  via  minor  modifications  to  these  GOFT.®. 

1.2  Problem  Statement 

The  Anderson-Darling  (A-D)  and  Cramer-von  Mises  (C-VM)  GOFT  statistics 
have  been  modified  by  earlier  researchers  to  be  more  powerful  for  the  normal  distri¬ 
bution;  however,  I  believe  that  there  is  room  for  improvement.  My  research  goal  is 
to  modify  the  A-D  and  C-VM  tests  to  produce  a  better  GOFT  for  normality  with 
the  parameters  estimated  from  the  sample. 

i 

l.S  Methodology  ! 

The  research  effort  will  consist  of  the  following  steps: 

1.  Produce  unmodified  A-D  and  C-VM  critical  value  tables  for  the  normal  distri- 

i 

bution. 

2.  Produce  modified  A-D  and  C-VM  critical  value  tables  for  the  normal  distribu¬ 
tion. 

3.  Conduct  tests  to  compare  the  powers  of  these  modified  and  unmodified  statis¬ 
tics. 


l.S.l  Critical  Value  Tables.  In  the  first  step,  critical  value  tables  will 
be  generated  for  the  unmodified  A-D  and  C-VM  statistics.  These  tables  will  be 
developed  by  generating  random  samples  from  a  normal  distribution  via  Monte  Carlo 
simulation  using  samples  of  size  four. 

Then  using  the  appropriate  equations  outlined  in  chapter  three,  the  A-D  and 
C-VM  statistics  will  be  calculated.  These  steps  will  be  repeated  10,000  times,  using 
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a  sample  size  of  four,  to  generate  10,000  independent  A-0  statistics  and  10,000 
independent  C-VM  statistics.  These  statistics  will  then  be  used  to  interpolate  the 
critical  values  for  the  respective  tables.  Once  this  is  accomplished,  it  will  then  be 
repeated  for  samples  of  size  n  =  5,6,... ,50. 

This  same  sequence  of  events  will  be  accomplished  using  the  modified  A-D  and 
C-VM  statistical  values.  Which  will  in  turn  allow  step  three  to  be  worked  as  we  will 
have  critical  value  tables  for  both  the  unmodified  and  the  modified  A-D  and  C-VM 
tests. 


/.5.J?  Power  Comparison.  The  power  of  a  statistical  test  is  the  probability 
of  correctly  rejecting  a  false  null  hypothesis.  The  gist  of  the  power  study  is  in  toting 
to  see  if  the  null  hypothesis  is  rejected  when  in  fact  it  should  be  rejected;  i.e.,  what 
percentage  of  the  time  is  the  null,  that  I  have  a  normal  sample,  rejected  when  in 
fact  my  sample  ^  from  an  exponential.  The  higher  the  percentage,  the  better,  more 
powerful  the  test.  These  powers  of  the  modified  and  unmodified  A-D  and  C-VM  tests 
will  be  compared  to  determine  which  amongst  them  can  best  detect  whether  or  not 
a  sample  of  data  comes  from  a  population  characterized  by  the  normal  distribution. 

To  perform  these  power  tests,  the  A-D  and  C-VM  test  statistics  will  be  cal¬ 
culated  under  the  null  hypothesis  that  the  random  samples  follow  the  normal  dis¬ 
tribution  with  the  specified  parameters.  The  calculated  statistics  for  each  test  will 
then  be  compared  to  the  appropriate  row  and  column  of  the  corresponding  critical 
value  tables  obtsuned  in  steps  one  and  step  two  to  determine  whether  to  reject  the 
null  hypothesis.  This  wili  be  repeated  10,000  times.  Then,  the  number  of  times  each 
statistic  exceeds  the  respective  critical  value  will  be  divided  by  10,000  and  this  wili 
give  us  the  power  of  the  test  -  the  percentage  of  the  time  the  null  hypothesis  is  cor¬ 
rectly  rejected.  This  power  computation  will  be  performed  with  data  created  from 
Monte  Carlo  simulation  for  the  uniform,  exponential,  double  exponential,  Weibull, 
beta  and  the  log-normal  distributions.  Upon  completion  of  power  computations,  the 
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powers  of  the  modified  A-D  and  C-VM  GOFTs  will  be  compared  to  the  unmodified 
A-D  and  C-VM  GOFTs  to  determine  which  test  is  most  powerful. 
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II.  Literature  Review. 


2.1  Introduction 

I  have  accomplished  a  literature  review  on  the  goodness-of-fit  test  techniques. 
Research  studies  have  been  investigated  through  the  1991  Current  Index  to  Statistics 
(the  most  current  available),  the  Defense  Technical  Information  Center  resources, 
and  Air  Force  Institute  of  Technology  theses.  My  discussion  on  hypothesis  testing, 
test  statistics,  critical  values,  the  bootstrap  technique,  plotting  positions,  power 
studies,  and  Monte  Carlo  simulation  is  important  because  these  things  are  crucial 
to  the  conclusion  of  this  thesis. 

2.2  Goodness-of-fit  Tests 

Statistics  is  much  more  art  than  science.  While  scientists  and  engineers  of¬ 
ten  devise  and  continually  use  statistical  methods  in  their  work,  statistics  is  still  a 
subjective  art  (17).  The  field  of  Statistics  has  many  useful  bramches  and  procedures 
within  it.  As  mentioned  earlier,  one  of  these  areats  hats  been  in  verifying  a  distribution 
for  a  set  of  observations.  Over  the  yeau's  mainy  different  test  procedures  have  been 
developed  to  au;complish  this  amd  they  have  come  to  be  known  as  goodness-  of-fit 
tests  (7).  The  data  are  assumed  to  be  saunples  from  the  hypothesized  distribution 
unless  there  is  sufficient  evidence  to  disprove  the  assumption  (1:72). 

Similau'  to  all  other  statistical  work,  there  exist  potentiad  problems  with  the 
goodness-of-fit  tests.  One  possible  problem  is  that,  even  though  one  test  shows  there 
is  insufiBcient  evidence  to  disprove  that  the  saunple  came  from  a  certaiin  distribution, 
there  may  be  other  distributions  which  fit  the  data  better  if  only  they  aire  tested. 
Or  it  may  not  be  am  entirely  different  distribution,  but  the  saune  distribution  with 
different  parauneters.  As  such,  it  is  importaint  to  remember  that  the  goodness-of-fit 
test  only  rejects  the  assumed  distribution  when  there  is  sufficient  evidence  that  the 
distribution  is  incorrect;  it  does  not  mean  that  the  selected  distribution  is  the  best 
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one  (1:73).  In  order  to  perform  a  goodness-of-fit  test,  it  is  important  to  ha^’e  an 
understanding  of  how  a  hypothesis  test  works. 


2.3  Hypothesis  Testing. 

Hypothesis  testing  is  the  process  of  inferring  from  a  sample  whether  to  accept 
a  certain  statement  about  the  population  from  which  the  sample  was  taken(5:76). 
Any  statistical  test  of  a  hypothesis  works  in  this  way  (16:429): 


1.  The  null  hypothesis,  Hq,  is  stated  in  terms  of  the  population. 

2.  The  alternative  hypothesis,  is  stated  in  terms  of  the  population. 

3.  A  test  statistic  is  chosen. 


4.  A  rule  is  made,  in  terms  of  possible  values  of  the  test  statistic,  for  deciding 
when  to  accept  the  null  hypothesis  and  when  to  reject  it  -  the  rejection  region. 

5.  Based  on  a  random  sample  from  the  population,  the  test  statistic  is  evaluated, 
and  a  decision  is  made  to  accept  or  reject  the  null  hypothesis. 


The  test  statistic  is  a  function  of  the  sample  data  upon  which  the  statistical 
decision  will  be  based.  If  for  a  particular  sample  the  computed  value  of  the  test 
statistic  fdls  in  the  rejection  region,  the  null  hypothesis  Hq  is  rejected  and  the 
alternative  hypothesis  Ha  is  aiccepted.  If  the  value  of  the  test  statistic  does  not  fall 
into  the  rejection  region,  then  Ho  is  ju:cepted  (16:429)(5:77-78).  Because  l^tatistics 
is  not  an  exact  science,  it  is  not  surprising  that  there  are  two  ways  of  rhaking  a 
mistake  in  a  hypothesis  test.  They  are,  either  rejecting  a  true  null  hypothesis,  or 
accepting  a  false  null  hypothesis.  These  two  error  types,  types  I  and  II  respectively, 
have  associated  with  them  certain  probabilities  of  their  being  made,  designated  as 
a  and  ^  (5:  79). 

The  following  example  illustrates  how  a  statistic2d  test  works. 
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A  psychological  study  was  conducted  to  compare  the  reaction  times 
of  men  and  women  to  a  certain  stimulus.  Independent  random  seimples 
of  50  men  and  50  women  were  employed  in  the  experiment.  The  results 
are  shown  in  table  2.1.  Do  the  data  present  sufficient  evidence  to  suggest 
a  difference  between  true  mean  reaction  times  for  men  and  women?  Use 
a  =  .05. 


Table  2.1.  Hypothesis  Test  Example  Data 


Men 

Women 

ni  =  50 

0 

11 

yi  =  3.6  seconds 

yi  =  3.4'}  seconds 

te 

II 

00 

II 

Solution:  Let  /ij  and  denote  the  true  mean  reaction  times  for  men 
and  women,  respectively.  Then  if  we  wish  to  test  the  hypothesis  that  the 
means  are  equal,  we  will  test  Ho  :  (fit  —  ^3)  =  0  against  Ha  :  (^1  —  fii)  ^ 
0.  Note  that  we  use  the  two-sided  alternative  to  detect  either  the  case 
Hi  <  Hi  or  the  reverse  Hi  >  io  case  Ho  is  false. 

The  point  estimator  of  {hi  ~~  Hi)  Js  it  satisfies  the 

assumptions  of  our  large-sample  test.  Hence  if  we  desire  to  test  Ho  : 
Hi~Hi  —  ^0  (Do  fixed)  versus  any  alternative,  the  test  statistic  is  given 
by 


—  I^)  —  Do 


ni  nj 


(2.1) 


where  aj  and  otj  are  the  respective  population  variances.  In  this 
application  we  desire  a  two-t£uIed  test.  Thus  for  a  =  .05  we  reject  Ho  for 
\z\  >  z^=  Z.02S  =  1.96 

For  large  samples  (say,  n  >  30)  the  sample  vari2mces  provide  good 
estimates  of  their  corresponding  population  variances.  Substituting  these 
values,  along  with  jfTjl/aj  and  Do  =  0,  into  the  formula  for  the  test 

statistic,  we  have  z  =  ^  =  —2.5 

VW  +  M 

This  value  is  less  than  —z^  =  —1.96  and  therefore  falls  in  the  rejection 
region.  Hence  we  reject  the  hypothesis  of  no  difference  in  mean  reaction 
times  for  men  and  women  (16:439-440). 


There  are  many  different  types  of  hypothesis  tests,  one  of  the  more  common 
being  of  the  type  of  the  above  example.  This  review  is  concerned  with  a  more  specific 
type  of  hypothesis  test  -  the  goodness-of-fit  test.  The  goodness-of-fit  tests  require  a 
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slightly  different  test  statistic  from  the  Z  in  the  example  above,  which  leads  us  to 
the  next  section  -  Test  Statistics. 

S.4  Test  Statistics 

Two  major  classes  of  test  statistics  used  for  goodness-of-fit  tests  are  the  Chi- 
square  and  the  Empirical  Distribution  Function  statistics. 

2.4.1  Chi-square  Test  Statistics.  The  first  goodness-of-fit  test  to  be  devel¬ 
oped  was  the  Chi-square  test  created  by  Karl  Pearson  in  1900.  After  working  with 
data  analysis  for  many  years,  he  gave  up  on  the  notion  that  biological  populations 
were  normally  distributed,  he  created  other  distributions,  and  then  created  his  Chi- 
square  test  to  determine  which  of  these  distributions  a  sample  came  from.  Modern 
developments  have  increased  the  flexibility  of  the  Chi-square  tests,  especially  when 
unknown  parameters  must  be  estimated  in  the  hypothesized  family  (7:63).  The 
motivating  force  of  the  Chi-square  test  is  that  it  compares  the  observed  frequen¬ 
cies  with  expected  frequencies  of  the  hypothesized  cumulative  distribution  function. 
This  test  is  restricted  to  large  samples,  approximately  25  or  greater  (16:310).  Even 
with  the  development  of  other  goodness-of-fit  tests,  the  advantages  of  the  Chi-square 
test  are  not  diminished;  it  is  an  economical  first  cut  at  trying  to  figure  the  samples 
distribution,  it  C2m  be  applied  to  discrete  populations,  and  it  can  be  modified  when 
parameter  values  are  unknown  by  reducing  the  number  of  degrees  of  freedom  (15:68). 

2.4-2  Empirical  Distribution  Function  (EDF)  Statistics.  In  addition  to  the 
basic  Chi-square,  there  is  the  general  class  of  statistics  called  Empirical  Distribution 
Function  statistics  which  are  £dso  used  for  the  goodness-of-fit  test.  The  empiricail 
distribution  function  is  a  step  function,  taking  a  value  between  0  and  1  for  each 
member  of  the  ordered  sample,  to  estimate  the  cumulative  distribution  function. 
EDF  statistics  are  just  like  their  Chi-square  counterpwts  in  that  they  measure  the 
discrepancy  between  the  EDF  and  the  hypothesized  cumulative  distribution  function 
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(7:  97).  For  a  given  random  sample  of  size  n,  let  be  the  ordered 

statistics.  Suppose  further  that  the  distribution  of  X  is  F(x).  Then,  the  empirical 
distribution  function  is  F„(x)  defined  by 

.  .  number  of  observations  <x  /«  o\ 

Fn{x)  = - - -  (2.2) 

where 

X  ranges  between  — oo  and  oo 

More  precisely,  the  definition  becomes  (7:98-99): 

■fn(®)  =  0,  X  <  A'(i)  (2.3) 

Fn{x)  =  X(i)  <  X  <  X(,+i),  *  =  1,  ...,n  -  1  (2.4) 

where  i  is  the  rank  of  Xi  amongst  all  the  ordered  observations. 

F„(X)  =  1,  ^(n)<X  (2.5) 

F„(x)  is  a  step  function  that  t£ikes  a  step  up  of  height  1/n  as  each  ordered 
sample  observation  is  reached.  One  type  of  EDF  statistic  is  the  Kolmogorov- Smirnov 
statistic. 

2.4.S  The  Kolmogorov-Smimov  (K-S)  Statistic.  The  K-S  statistic  is  an 
EDF  statistic  in  wide  use.  It  is  the  largest  vertical  distance  between  the  completely 
specified  hypothesized  cumulative  distribution  function  (CDF)  and  the  observed  em¬ 
pirical  distribution  of  the  sample  data(7:204).  When  it  is  necessary  to  completely 
specify  all  the  parameters  of  the  distribution,  then  use  of  the  K-S  statistic  is  appro¬ 
priate,  but  only  if  the  entire  specification  can  be  made  prior  to  sampling  (3:296). 
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Work  done  by  D’Agostino  and  Stephens  has  shown  that  EDF  statistics  are  much 
more  powerful  thcin  Pearson  Chi-  square  statistics  (7:110).  Another  EDF  statistic 
mentioned  in  the  previous  chapter  is  the  C-VM  statistic. 

2.4.4  The  Cvamer-von  Misea  Statistic.  The  C-VM  statistic  is  based  on  the 
squared  summation  of  the  difference  between  the  EDF  and  the  distribution  being 
tested.  This  test  statistic  is  defined  as  (7:101): 

=  (2.6) 

where  n  =  sample  size, 

Uf  =  F(Ar(,))  =  CDF  for  the  distribution  of  interest, 
i  is  the  rank  of  AT,-  amongst  all  the  ordered  observations, 
aod  -^(i)  ^  Af(2)  <  ...  <  are  the  n  observations  in  order. 

The  C-VM  test  statistic,  like  the  K-S  test  statistic,  is  also  a  function  of  the 
vertical  distance  between  the  hypothesized  distribution  and  the  Empiriczd  Distribu¬ 
tion  Functions.  They  differ  in  the  fact  that  the  C-VM  test  goes  one  step  further  by 
considering  every  difference  between  the  curve  formed  by  the  CDF  and  the  curve 
formed  by  the  EDF.  Conover  believes  that  the  C-VM  test  is  more  appealing  than 
the  X-S  test  because  the  C-VM  statistic  uses  more  of  the  sample  data.  (5:306). 

The  following  is  a  numerical  example  of  how  a  C-VM  goodness-of-fit  test  is 
performed  on  a  sample  of  data,  n=8: 

Ho‘  The  sample  data  in  table  2.2  came  from  a  normal  distribution 
with  fi  =  4.3  and  <7  =  3.039. 

Ha’.  The  sample  came  from  some  other  population. 

Decision  Rule:  If  the  C-VM  statistic  >  .1182  reject  the  null  hypoth¬ 
esis,  otherwise  2K:cept  the  null  hypothesis,  (note:  using  a  =  .05  and  n=8 
the  C-VM  critical  value  from  the  table  on  page  4-3  :3  .1182) 
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The  cumulative  distribution  function  for  the  normal  distribution  be¬ 
ing: 

f (x)  =  f  f(x)Jx  /■  _•  (2.7) 

J-oo  J-OO  ay/2K 

And  remember  that  C  -  VAf  =  +  Er=i(£^.  - 


Table  2.2.  C-VM  GOFT  Example  Data 


j 

X 

F(x)=G, 

A  -  ^ 

-  2n 

(Vi  -  A,y 

1 

0.2 

.1038 

.0625 

.00171 

2 

1.6 

.2033 

.1875 

.00025 

3 

2.1 

.2483 

.3125 

.00412 

4 

3.0 

.3446 

.4375 

.03863 

5 

4.8 

.5596 

.5625 

.00001 

6 

5.0 

.5871 

.6875 

.01008 

7 

8.1 

.8790 

.8125 

.00442 

8 

9.6 

.9484 

.9375 

.00012 

G  -  VM  =  ^  +  .02934  =  .01042  +  .02934  =  .03976 
Because  .03976  2  .1182,  we  accept  the  null  hypothesis  that  the  sample 
does  indeed  come  from  a  normal  distribution  with  ft  =  4.3  and  <r  =  3.039 
(22:17-18). 

I 

Another  EDF  statistic  mentioned  in  the  previous  chapter  is  the  A*D  statistic. 

! 

2.4,5  The  Anderson- Darling  Statistic.  Anderson  and  Darling  created  a 

I 

new  statistic  by  incorporating  a  weight  function  into  the  K-S  and  C-VM  statistics. 
The  A-D  test  statistic  is  defined  as  (7:100-101): 


A-D  =  -n  - 2f; -  K..,,,))  (2.8) 

i=t  2n 

where  n  =  sample  size, 

Ui  =  F(A’(,))  =  CDF  for  the  distribution  of  interest, 
i  is  the  rank  of  A’,-  amongst  all  the  ordered  observations, 
and  X^l)  <  A(2)  <  ...  <  X^„)  are  n  observations  in  order. 
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The  A-D  statistic  is  based  on  a  weighted  average  of  the  squared  discrepancies 
between  the  curve  formed  by  the  CDF  auid  the  curve  formed  by  the  EDF. 

The  following  is  a  numerical  example  of  how  an  A-D  goodness-of-llt  test  is 
performed  on  a  sample  of  data,  n=8: 

Hq'.  The  sample  data  in  table  2.3  came  from  a  normal  distribution 
with  n  =  4.3  and  a  —  3.039. 

Ha,',  The  sample  came  from  some  other  population. 

Decision  Rule:  If  the  A-D  statistic  >  .6549  reject  the  null  hypothesis, 
otherwise  accept  the  null  hypothesis,  (note:  using  a  =  .05  and  n=8  the 
A-D  critical  value  from  the  table  on  page  4-2  is  .6549) 

Remember  that  A-D  =  -n-  2  ^(loge(<^*)  +  ^05e(l-tAt-i+i)) 

and  F(x)  =  /f<„/(x)dx  =  Also,  Ai  =  ln(f/.)  and 

Hi  =  ln(l  - 

Table  2.3.  A-D  GOFT  Example  Data 


rr 

X 

F(x)=t;,- 

Un-i-t-l 

(2i-l)(Ai-fBi) 

1 

0.2 

.1038 

.9484 

-2.265 

-2.964 

-5.229 

2 

1.6 

.2033 

.8790 

-1.593 

-2.112 

-11.115 

3 

2.1 

.2483 

.5871 

-1.393 

-.885 

-11.390 

4 

3.0 

.3446 

.5596 

-1.065 

-.820 

-13.195 

5 

4.8 

.5596 

.3446 

-.581 

-.423 

-9.036 

6 

5.0 

.5871 

.2483 

-.533 

-.285 

-8.998 

7 

8.1 

.8790 

.2033 

-.129 

-.227 

-4.628 

8 

9.6 

.9484 

.1038 

-.053 

-.110 

-2.445 

EJ=.i(2;-1)(A^-|-D^)  =  .02934 

A  -  D  =  -8  -  (|)(-66.036)  =  -8  -f-  8.2545  =  .2545  Because  .2545 
.6549,  we  accept  the  null  hypothesis  that  the  saunple  does  indeed  come 
from  a  normal  distribution  with  /x  =  4.3  and  c  =  3.039  (22:16-17). 

The  critical  values  used  in  the  two  tests  above  were  taken  from  critical  value 
tables,  which  leads  to  the  next  pwt  of  this  review,  critical  values. 


2.5  Critical  Values. 

Previous  critical  values  have  been  generated  by  taking  random  samples  from 
the  test  distribution,  calculating  the  test  statistic  and  then  ordering  the  statistics 


2-8 


from  smallest  to  largest.  The  best  way  to  explain  the  nature  of  a  critical  value  is 
with  a  word  picture:  if  100  independent  r<indom  samples  of  size  n  are  generated,  and 
from  them  100  goodness-of-fit  statistics  are  calculated  and  ordered  from  smallest 
to  largest,  then  the  critical  value  for  1  —  or  =  1  —  .05  =  .95  is  the  95‘*  largest 
order  statistic.  This  reasoning  works  well  when  there  is  a  one-to-one  correspondence 
between  the  a-level  and  the  order  statistic,  but  what  if  for  the  example  with  100 
statistics  a  critical  value  was  needed  at  a  =  .0005?  There  would  be  a  problem  because 
in  the  example,  discrete  values  are  being  used  to  determine  critical  values  for  what 
is  actually  a  continuous  distribution.  To  get  around  this  conflict,  all  that  must  be 
done  is  an  interpolation  of  the  data,  and  for  endpoints,  an  extrapolation.  This  is 
simply  a  method  for  representing  these  order  statistics  on  a  continuous  spectrum. 
This  is  done  by  plotting  the  values  of  the  order  statistics  and  representing  the  spaces 
between  them  as  piecewise  linear  functions.  To  6nd  a  critical  value,  all  that  is 
necessary,  graphically,  is  to  find  the  value  on  the  x-axis  that  corresponds  with  the 
selected  1  —  o  value  of  Interest  on  the  y-axis  (8). 

First,  the  n  ordered  test  statistics,  11,0:3,  are  plotted  against  their  1— a  = 
yitVit  •••?  Vn  which  are  calculated  with  the  median  rank  formula  of  the  form:  y,-  = 
i=the  rank  of  the  ordered  statistic  and  n=s£imple  size.  Letting  the  ordered  statistics 
be  represented  by  the  horizontal  axis  and  letting  the  y,-  be  represented  by  the  vertical 
axis  and  not  coincidentally  be  scaled  between  zero  and  one,  the  statistics  will  be 
represented  on  a  continuous  function.  After  this  is  accomplished,  the  endpoints 
are  found  by  linear  extrapolation.  The  first  point  on  the  horizontal  axis,  xq,  is 
computed  by  linearly  extrapolating  from  the  second  and  third  points,  subject  to  a 
non-negativity  restriction.  Extrapolation  is  performed  by  using  the  standard  linear 
slope-intercept  formula  y  =  m*x-\-bto  compute  the  endpoints  xo  and  x„+i.  To  find 
the  first  endpoint  on  the  horizontal  axis,  the  slope  is  calculated  by: 

_ Vi+i  -  Vi 

m  ~ - 

®«+i  ~ 


(2.9) 
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and  the  intercept  is  found  by: 


b  =  yi  —  m*  Xi  (2-10) 

Then  the  lower  endpoint  iq  is  found  by: 

X  =  _  Q  ~  ^ 

m  mm 

Using  the  non-negativity  rule, 

xo  =  >7iax(0, — -)  (2.12) 

m 

After  both  endpoints  are  calculated,  n  test  statistic  values  and  their  two  ex¬ 
trapolated  endpoints  are  plotted  on  the  x  -  y  axis.  Then,  the  critical  values  corre¬ 
sponding  to  the  desired  percentiles  are  found  by  linearly  interpolating  between  the 
points  (xi,y,)  and  (x,+i,y,+i)  using  the  formulas  above;  where  Cp  is  the  critical  value 
for  the  desired  percentile  (20). 

To  find  a  criticcil  value,  all  that  is  necessary,  graphically,  is  to  find 
1  —  Q  on  the  vertical  axis  and  extend  along  the  line,  ?(,•)  =  1  —  a,  to 
intercept  the  plotted  function.  The  value  of  the  horizontal  component  is 
the  critical  value  of  the  statistic  at  significance  level  or. 

Finding  the  critical  value  with  a  computer  requires  finding  the  largest 
Y{i)  that  is  less  than  1— q.  Suppose  that  Y(k)  is  the  largest  rank.  Then, 
the  standard  lineu  slope-intercept  formula  {y  —  mx  -|-  h)  is  used  to  find 
the  critical  value.  The  chjuige  in  y  can  be  found  using  Y(k)  and  Y^k+i). 
Similarly,  the  change  in  x  can  be  found  using  X{^k)  and  ^(*+1).  After 
finding  the  con8t2Lnt,  b,  at  (A’(*),y(i)),  one  can  then  let  y  equal  1  —  q  in 
order  to  find  x,  the  critical  value. 

Suppose  ten  samples  are  taken  as  shown  in  table  2.4: 

In  figure  2.1,  the  statistics  aure  plotted  versus  their  median  step  ranks. 

And  from  table  2.4,  Y^lf  =  .05  7(2)  =  -15  Xm  =  .22  ^(2)  =  .41  Using 
the  equation,  y  =  mx  -f  6,  m  =  slope  =  =  .505  and 

b  =  Y(i)  -  mX^l)  -  =  .067  -  (.505)(.22)  =  -^.0441  and  i  =  ^  =  ^  = 


(2.11) 
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Table  2.4.  Critical-value  Example  Data 


i 

Median  Rank  yi  = 

Statistics  X,- 

1 

.067 

.22 

2 

.163 

.41 

3 

.26 

.42 

4 

.356 

.67 

5 

.452 

.98 

6 

.548 

1.02 

7 

.644 

1.03 

8 

.74 

1.08 

9 

.837 

1.12 

10 

.933 

1.13 

m 


Figure  2.1.  Critical  Value  Example 


f 


.0873  Since  x  =  .0873  >  0,  A'(o)  =  x.  Again,  if  x  had  been  less  than  zero, 

^(0)  would  have  been  set  to  zero.  Extrapolation  for  X(n)  is  performed 
the  same  way.  V/jo),  V(9),  X(io),  and  JV(9)  are  used  to  find  the  slope.  The 
constant,  b,  is  calculated  at  either  (X(tc),  VJio))  or  (X(9),y(9)).  Then 
^(11)  =  where  m  is  the  slope.  Now  that  the  function  is  continuous 
(by  extrapolation)  on  the  intervad  (0,1),  the  critical  values  can  be  found. 

At  o  =  .10,  previous  studies  would  have  picked  1.12,  or  the  9‘^  largest 
statistic  as  the  critical  value.  Using  the  bootstrap  method,  the  value  is 
1.1266  (if  the  given  median  step  rank  is  used).  To  get  the  critical  value 
using  the  bootstrap  technique,  the  largest  TJ,)  less  than  or  equal  to  .90 
is  found.  In  this  case,  this  is  yjg)  =  .85.  Therefore,  k=9  and  +  1  =  10. 
Tk'".  "■  =  =  rSem  =  ^  =  9-6  6  =  .837  -  (9,61(1.12)  = 

.837  -  10.752  =  -9.915  and  crilicol  value  =  =  1.1266  As 

VaO 

one  can  see,  the  critical  value  will  vary  with  statistics  calculated  from 
random  samples.  (22:19-23). 

Typically,  researchers  are  interested  in  getting  the  critical  values  for  several 
o-levels  such  as  a  =  .01, ...,  .20.  Once  these  critical  ’^alue  tables  have  been  generated 
and  tabled  for  a  particular  distribution,  an  assessment  of  the  usefulness  of  the  table 
needs  to  be  performed.  This  assessment  has  come  to  be  known  as  a  power  study. 


2.6  Power  Studies. 

The  power  of  the  goodness-of-fit  test  is  the  probability  that  the  test  will  reject 
the  null  hypothesis  when  in  fact  it  should  reject  the  null  hypothesis  (16:473).  The 
statistician  wants  his  test  to  be  as  powerful  as  possible  for  identifying  those  estimates 
of  distributions  that  are  flat-out  wrong.  The  following  discussion  presents  some  of 
the  work  done  in  efforts  to  create  more  powerful  goodness-of-fit  tests. 

Green  and  Hegazy  teamed-up  to  attempt  to  improve  the  standard  EDF  goodness- 
of-fit  tests.  They  modified  the  A-D,  C-VM,  and  the  K-S  goodness-of-fit  tests  where 
the  parameters  were  estimated  from  the  sample.  They  showed  that  the  modified 
A-D  and  C-VM  tests  achieved  higher  power  (9:204). 

Porter  modified  the  A-D,  C-VM  and  the  K-S  goodness-of-fit  tests  for  the  three- 
parameter  Pareto  distribution.  He  concluded  that  the  power  of  his  modified  tests 
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improved  as  the  sample  sizes  increased.  He  also  showed  that  A-D,  C-VM  and  K-S 
tests  were  more  powerful  than  the  Chi-square,  especially  when  the  sample  data  were 
taken  from  the  Weibull,  beta,  or  norma]  distributions  (20:7). 

Ream  modified  the  K-S,  A-D,  and  C-VM  tests  for  the  normal  distribution  to 
test  the  technique  of  reflecting  data  points  about  the  mean  iuid  created  tables  of 
critical  values  of  the  modified  K-S,  A-D,  and  C-VM  statistics.  He  concluded  that 
against  symmetric  distributions  the  power  increased,  whereas  against  non-symmetric 
distributions  a  significant  loss  of  pov/er  occurred.  He  further  recommended  that  the 
consequences  of  this  fact  be  strongly  considered  before  use  of  his  new  technique  and 
statistics  (22:62-63). 

Bush  modified  the  A-D  and  C-VM  tests  for  the  Weibull  distribution.  He  pre¬ 
sented  a  Monie  Carlo  method  for  obtaining  the  critical  values  of  the  modified  A-D 
and  C-VM  goodness-of-fit  tests  for  the  three  parameter  Weibull  distribution  when 
the  scale  and  location  parameters  were  not  specified.  He  also  concluded  that  A-D 
and  C-VM  tests  were  not  very  powerful  when  the  sample  size  was  five.  When  the 
sample  size  increased,  he  showed  that  A-  D  and  C-VM  were  more  powerful  than  K-S 
tests  (4:52). 

Viviano  modified  the  A-D  and  C-VM  tests  for  the  Gamma  distribution.  He 
showed  that  the  A-D  goodness-of-fit  test  was  more  powerful  against  the  log-normal 
distribution  than  the  C-VM  test(24:47). 

Kahya  modified  the  A-D  and  C-VM  tests  by  adjusting  the  plotting  positions  of 
the  original  statistics.  He  showed  that  the  power  of  the  A-D  test  increased  for  most 
all  distributions  and  the  power  of  the  C-VM  test  increased  for  the  uniform,  Weibull, 
exponential,  double  exponential,  and  log-normal  (13:6-1). 

According  to  Stephens,  critical  vadue  tables  generated  from  one  continuous 
distribution  can  be  applied  to  GOFTs  for  all  continuous  distributions  as  long  as  the 
parameters  are  completely  specified  for  the  sample  (7:102). 


'  i  /  '  .  '  ■  ■/ 
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All  work  with  GOFTs  involves  many  tests.  In  fact,  to  do  work  of  any  signifi¬ 
cance,  many  tests  must  be  performed  on  samples  of  data.  And  that  leads  us  to  the 
last  topic  of  this  chapter,  Monte  Carlo  simulation. 


2.1  Monte  Carlo  Simulation 

In  many  cases  a  system  or  phenomenon  under  study  is  either  too  complex  to 
be  adequately  defined  by  mathematical  formulation  or  the  mathematicaJ  formula¬ 
tion  is  extremely  difficult,  or  even  impossible,  to  solve  by  analytic  techniques.  When 
these  difficulties  occur,  a  solution  can  often  be  obtained  by  use  of  the  technique 
mentioned  earlier  -  simulation.  Simulation  is  a  procedure  for  imitating  (not  repre¬ 
senting)  a  system  or  phenomenon  for  different  values  of  the  controllable  variables. 
Much  of  the  work  with  GOFTs  use  Monte  Carlo  simulation  to  create  data  that 
mimics  many  different  populations:  e.g.  normal,  log-normal,  exponential,  double 
exponential,  uniform,  and  Weibull.  It  is  one  thing  to  have  data  that  is  representa¬ 
tive  of  a  distribution,  e.g.  ten  samples  from  an  exponential  distribution  with  A  =  10, 
but  it  is  quite  amother  to  have  that  sample  be  randomly  generated  again  and  again 
and  again  (17).  The  approach  is  to  observe  random  variates,  chosen  so  that  they 
directly  simulate  the  random  processes  of  the  original  problem.  Then  based  upon 
the  multiple  observations,  the  desired  general  solution  is  inferred  from  the  behavior 
of  the  random  numbers  (10:2-4).  The  main  weakness  in  Monte  Carlo  simulation  is 
that  the  answers  it  produces  are  to  some  degree  uncertadn  since  they  are  inferred  raw 
observational  data  consisting  of  raw  numbers.  This  weakness  must  be  accounted  for 
because: 

Whenever  one  is  inferring  general  laws  on  the  basis  of  particular  obser¬ 
vations  associated  with  them,  the  conclusions  are  uncertain  inasmuch  as 
the  particular  observations  are  only  a  more  or  less  representative  sample 
from  the  totality  of  all  observations  which  might  have  been  made.  Good 
experimentation  tries  to  ensure  that  the  sample  shall  be  more  rather  than 
less  representative...  (Monte  Carlo  answers]  can  nevertheless  serve  a  use¬ 
ful  purpose  if  we  can  manage  to  make  the  uncertainty  fairly  negligible. 
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that  is  to  say  to  make  it  unlikely  that  the  answers  are  wrong  by  very 
much  (10:4-5). 

Thus  there  is  usually  no  cause  for  concern  if  the  uncertainty  is  negligible  for 
practical  purposes.  One  way  of  reducing  uncertainty  is  to  base  the  Monte  Carlo  work 
on  a  larger  number  of  observations.  This  is  done  because  the  law  of  large  numbers 
states  that,  as  the  sample  size  increases,  the  differences  between  the  sample  mean 
and  the  population  mean  becomes  smaller  (1:176). 

All  work  done  in  this  area  has  not  been  done  to  satisfy  one  particular  seimple 
of  data  from  say  the  V/eibull  distribution.  But  instead,  the  researchers  were  trying 
to  establish  general  principles  that  applied,  if  not  across  all  probability  distribution 
families,  then  at  least  across  one  distribution  family  -  regardless  of  sample  size.  To 
achieve  this  end,  Monte  Carlo  simulation  was  key. 

2.8  Conclusion. 

Every  single  one  of  these  research  studies  has  been  done  in  order  to  try  to 
find  more  powerful  tests  for  the  different  distributions.  They  all  modified  the  test 
techniques  or  used  a  new  technique  for  only  one  purpose;  to  obtain  a  "better”  test 
so  that  the  analysts  are  better  able  to  determine  how  well  a  set  of  data  fits  some 
distribution.  However,  there  are  still  many  different  avenues  that  can  be  pursued 
to  find  even  more  powerful  tests  for  determining  goodness-of-fit.  Research  must  be 
done  in  this  area. 
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III.  Methodology 

5.1  Introduction 

This  chapter  describes  the  procedures  used  to  accomplish  the  research  objec¬ 
tives  of  this  thesis.  The  discussion  covers  the  modification  of  the  A-D  and  C-VM 
statistics,  calculation  of  the  critical  values,  and  finally,  calculation  of  the  power  ta¬ 
bles. 

5.2  Discussion 

The  research  effort  will  consist  of  the  following  parts: 

1.  The  A-D  and  C-VM  critic^d  value  tables  will  be  derived  using  the  known  A-D 
and  C-VM  statistic  formulas. 

2.  The  A-D  and  C-VM  critical  value  tables  will  be  generated  using  the  modified 
A-D  and  C-VM  statistic  formulas. 

3.  The  powers  of  the  modified  A-D  and  C-VM  statistics  will  be  determined  by 
Monte  Carlo  simulation  testing  against  alternative  distributions.  These  al¬ 
ternative  distributions  will  be  the  uniform,  exponential,  double  exponential, 
Weibull,  beta,  and  log-normal  distributions. 

4.  Analysis  will  be  conducted  on  the  powers  of  each  alternative. 

S.2.1  Modified  Test  Statistics.  As  stated  earlier,  the  work  of  this  thesis  is 
to  determine  which  A-D  and  C-VM  goodness-of-fit  statistics  are  most  powerful.  In 
addition  to  ex2unining  the  power  of  the  unmodified  A-D  and  C-VM  statistics,  this 
thesis  will  examine  five  modifications  of  each  of  the  A-D  and  C-VM  statistics  (a  total 
of  ten  new  statistics).  This  research  will  substitute  in  place  of  the  portion  of 
the  statistics  this  median  rank  quotient  Pi  =  ,  i=rank  of  the  ordered  statistic, 

and  n=8ample  size.  The  research  will  use  c  =  0,  0.3,  0.3175,  0.375  and  1  to  generate 
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these  five  different  median  rank  plotting  positions  for  development  of  the  modified 
A-D  and  C-VM  statistics: 


(3.1) 

(-0.3 
•  n  -f  0.4 

(3.2) 

,  t  -  0.3175 
*  "  n  -1-  0.365 

(3.3) 

,  t- 0.375 
*  n  -|-  0.25 

(3.4) 

n-1 

(3.S) 

Using  these  plotting  positions,  the  modified  A-D  and  C-VM  statistics  are: 


ADmedian,  c=o  =  — «  —  2  '  ,  1  Ooge(U.)  -I-  hge(l  —  Un_,+i))  (3.6) 

<=i  ”  ^ 


"  i  —  0  3 

ADmedian, 0=0.3  =  -n  -  2 ^  ^  ^  ^ (log,(f/,)  -H  logdl  -  Un-i+i)) 


ADmedian,  <7=0.3175 


”+ 


i  -  0.3175 
0.365 


(loge(U,)  -f  /o</e(l  -  Un-i+l  ))  (3-8) 
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ADmedun,c=o.37s  =  ^  -  Un-i+i))  (3.9) 

ADMEDtAKC=i  =  -n  ~  2]^  - — r(logj(t/,)  +  logdl  -  Un-i+i))  (3.10) 

1=1  “  * 

CVMmedian,c^o  =  4  (3-11) 

cvmmedian,  c=o.3  =  ^  Z(f'.  -  (3-12) 

CVA/M£DMAr.c=o.3,r5  =  ^  (3-13) 

CVMmedian,c=0273  =  ^  E(f^.  - 

CVMmedian,c=i  =  4  ^  ^(C^i  - 

S.S.2  Calculation  of  the  Critical  Values.  The  critical  Vcilue  tables  will  be 
generated  using  Monte  Carlo  simulation  of  random  deviates.  The  following  steps 
describe  the  procedure  for  generating  the  critical  value  table  for  the  unmodified  A-D 
when  the  sample  size  is  four. 

1.  Using  the  InternationaJ  Mathematical  and  Statistics  Library  (IMSL)  subrou¬ 
tine  RNNOR,  the  random  variates  will  be  generated  from  the  normeil  distribu¬ 
tion  with  a  mean  of  p  =  0  and  a  standard  deviation  of  it  =  1. 
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2.  The  random  variates  will  be  standardized  by  using  the  following  transforma¬ 
tion: 

=  fiZi  (3.16) 

3 

where  "S  =  the  sample  mean,  and  s  =  unbiased  sample 

standard  deviation. 

These  standau'dized  data  will  be  used  in  steps  three  through  six. 

3.  The  hjrpothesized  cumulative  distribution  function  for  the  normal  distribution, 
Ui,  will  then  be  calculated  for  i  =  1, 2, ...,  n. 

U,  =  /’  Hx)ix  =  /'  (3,17) 

J—oo  J—oo  3y2']f 


4.  The  A-D  statistic  will  be  calculated. 

5.  Steps  one  through  four  will  be  repeated  9,999  more  times  to  generate  a  total 
of  10,000  independent  A-D  statistics. 

6.  These  10,000  statistics  will  be  ordered  from  smallest  to  largest.  Using  the 
bootstrap  technique,  the  80‘*,85**,90‘^,95‘^,  and  99‘*  percentiles  of  the  distri¬ 
butions  of  the  A-D  statistic  will  be  determined  by  linear  interpolations.  These 
percentiles  represent  the  .20,  .15,  .10,  .05,  and  .01  levels  of  significance.  Then 
these  values  will  be  entered  into  the  A-D  critical  value  table. 

These  six  steps  will  be  done  again  for  n  =  5,...,  50  for  the  unmodified  A-D 
statistic.  Then  the  whole  procedure  will  be  repeated  for  the  five  modifications  to  the 
A-D  statistic,  the  unmodified  C-VM  statistic,  and  the  five  modified  C-VM  statistics 
as  well.  The  flowchart  representation  for  generating  the  critic2d  values  is  shown  in 
figure  3.1. 

S.2.S  Power  of  the  Goodness-of-fit  Tests.  Basicjdly,  the  power  of  a  goodness- 
of-fit  test  is  the  probability  that  the  test  will  reject  the  null  hypothesis  when  in  fact 
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Figure  3.1.  Generation  of  Critical  Values 
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it  should  be  reject^.  The  reason  for  conducting  the  power  test  is  to  help  determine 
if  the  modified  A-D  and  C-VM  statistics  result  in  more  powerful  goodness-of-fit  tests 
than  the  unmodified  A-D  and  C-VM  statistics. 


In  order  to  run  a  power  test,  random  deviates  from  non-Normal  distributions  of 
sample  size  n  will  be  generated  and  x  and  s  will  be  calculated.  Then,  the  A-D  and  C- 
VM  statistics  will  be  calculated  using  the  normal  CDF.  Because  of  the  earlier  work, 
the  critical  value  tables  are  available  for  look-up  of  critical  values  by  sample  size,  n, 
and  by  probability  of  a  type  I  error,  a.  The  calculated  A-D  and  C-VM  statistics 
will  then  be  compared  to  their  respective  table  value  to  see  if  the  test  passes  or  fails. 
This  process  will  be  repeated  10,000  times  and  the  number  of  times  each  statistic 
value  exceeds  the  respective  critical  v^Jue  will  be  counted  for  each  sample  size  of  n. 
This  total  number  will  represent  the  number  of  the  rejections  of  the  null  hypothesis. 

I 

Ideally,  this  count  would  be  10,000,  meaning  none  of  the  samples  passed  the  test  for 
coming  from  a  normsJ  distribution,  but  this  will  not  be  the  case  as  the  test  is  far 
from  clairvoyant.  This  total  rejection  number  will  be  divided  by  10,000  to  obtain 
the  power' 


The  power  = 


Total  number  of  rejections. 

10,000 


The  power  study  is  performed  at  n=10,  n=20,  n=30,  n=40,  and  n=50.After 
getting  all  the  powers  at  each  a-level,  final  conclusions  will  be  derived  for  e2w:h 
alternative  distribution. 


To  reiterate,  the  power  study  will  be  conducted  using  the  following  test  and 
the  steps  in  figure  3.2: 


Ho  :  Sample  deviates  follow  a  normal  distribution. 


Ha  :  They  follow  some  other  distribution. 
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Decision  Rule:  If  the  A-D  statistic  >  critical  value  reject  uhe  null  hypoth-^is, 
otherwise  accept  the  null  hypothesis 

S.S  Conclusion 

This  chapter  outlined  the  specific  steps  that  will  be  taken  to  generate  the  tests 
and  data  needed  to  accomplish  this  thesis.  The  next  chapter  will  show  the  tables 
generated  by  pursuing  the  methods  of  this  chapter. 
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Figure  3.2.  Generation  of  Power  tables. 
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IV.  Results 


Following  the  steps  outlined  in  the  previous  chapter,  this  chapter  presents  the 
critical  value  tables  and  the  power  study  tables. 

4.1  Critical  Value  Tables 

To  use  these  tables,  select  the  a  level  from  across  the  top  row,  and  then  go 
down  the  table  to  the  appropriate  sample  size  to  get  the  correct  criticcJ  value  for  the 
test. 
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Table  4.1.  Critical  VaJues  for  the  ANDERSON-DARLING  Test  with  Unmodified 
Plotting  Position;  alpha  =  .20  ...  .01 
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Critical  Values  for  the  CRAMER-VON  MISES  Test  with  Unmodified 
Plotting  Pwition;  alpha  =  .20  ...  .01 
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Table  4.3.  Critical  Values  for  the  ANDERSON-DARLING  Test  with  C  =  0.3175 
Median  Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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‘  T  Table  4.4.  Critical  Values  for  the  CRAMER- VON  MISES  Test  with  C  =  0.3175 

\  '  Median  Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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Table  4.5. 


Critical  Values  for  the  ANDERSON-DARLING  Test  with  C  =  0  Median 
Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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Table  4.6. 


Critical  Values  for  the  CRAMER-VON  MISES  Test  with  C  =  0  Median 
Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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Table  4.7.  Critical  Values  for  the  ANDERSON-DARLING  Test  with  C  =  0.375 
Median  Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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Critical  Values  for  the  CRAMER- VON  MISES  Test  with  C  =  0.375 
Median  Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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Critical  Values  for  the  ANDERSON-DARLING  Test  with  C  =  0.3  (Stan¬ 
dard)  Median  Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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Table  4.10. 


Critical  Values  for  the  CRAMER-VON  MISES  Test  with  C  =  0.3  (Stan¬ 
dard)  Median  Rank  Plotting  Position;  alpha  =  .20  ...  .01 


|K1 

imem 

.15 

1  -10 

1  .05 

.01 

In 

0.0759 

iBKsn 

0.1257 

rn 

IEESZ] 

IKiTTOyi 

I*.**!:!.] 

1  0.1031 

0.1403 

|n 

[eesb 

IKiTiTOTl 

Il2;ahl 

lElIES 

0.1500 

In 

IESSE] 

IEZ2Ea 

IEXS5S3 

0.1111 

lESmi 

rx 

igrron 

lES^ 

IliMU 

0.1133 

0.1544 

lEISSIl 

IMitUil 

lEZuUl 

0.1552 

11212^ 

iKssa 

1121^ 

0.1146 

0.1628 

lESSa 

lES^ 

0.1161 

KBsa 

lEZISIl 

lES^ 

112202 

iiAiaya 

IliSMI 

IE2^[ 

IQIQl 

msmM 

||2!Ea 

I21i2^ 

lEZiBia 

IliXlSil 

IESIE21 

lAELU 

■sisna 

iKiumi 

■ 

a 

■ 

ESSil 

IgilOCT 

1120^ 

iKiiigg 

IliZSiEi 

0.0995 

0.1211 

0.1714 

1  0.0874 

0.0995 

0.1221 

0.1770 

lllijiil 

0.0997 

0.1676 

\mma 

IKSIiSa 

0.1222 

0.1663 

EESiai 

IIBI 

IIKI 

[giltyigl 

0.1646 

t'lt'U-ia 

imi 

1130201 

0.1217 

0.1721 

■qniTiai 

mregti 

giinroii 

WiIgAM 

12!i£Qi 

ESEEal 

giign 

■■lliT3:7a 

■sraii 

■ntxirini 

■■HWgl 

0.1715  1 

msemm 

■iBUca 

Bsnii 

E£S1Z1 

■ifiT!?:?! 

0.1749 

EESll 

0.1241 

0.1756 

0.1217 

■tiTTO:;! 

mwiTiy 

0.1219 

0.1743 

0.1012 

0.1234 

0.1685 

■TTililiy 

esmi 

0.1225 

0.1701 

iiE2ia 

ESEZa 

0.1745 

■i8ij:W 

ESE2E1 

0.1237 

ESIQll 

liZMiJ 

mti!^ 

EIEZDi 

glUTWil 

Esssal 

■nfiCTra 

Ezumi 

Esm 

g»pynil 

EBEal 

EESilll 

mil 

ESEm 

ffnyyyi 

ESMlI 

Eisiiai 

I2EU1 

0.1231 

0.1779 

Essai 

EB3E1 

0.1227 

0.1728 

iiniigaii 

bftUiil 

inEa 

0.1787 

0.0788 

liAUUI 

I2ESI 

0.1741 

0.0797 

l2Ei:21 

0.1252  1 

0.1764 

ES:i3l 

raiijai 

rawai 

0.1727 

ni!^l 

raiiRiii 

LUiLII 

0.1736 

EEsai 

rarn™ 

MI’ISII 

ESESal 

ES^I 

L2mii 

0.1763 

EMEII 

usm9\ 

0.1231 

0.1772 

L2EUJI 

0.1253 

0.1734 

Table  4.11. 


Critical  Values  for  the  ANDERSON-DARLING  Test  with  C  =  1  Me¬ 
dian  Rank  Plotting  Position;  alpha  =  .20  ...  .01 
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4-2  Power  Study  Tables 

To  use  these  tables,  select  the  a  level  from  across  the  top  row,  and  go  down  the 
table  to  the  cluster  of  sample  sizes  of  interest  to  compare  the  six  different  powers. 


Table  4.13.  ANDERSON-DARLING  GOFT‘s  Power  Against  UNIFORM(0,1) 
data;  alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0,285 

0.220 

im 

0.079 

0.016 

Unmodified 

10 

0.286 

lEpi 

0.161 

0.081 

0.016 

C  =  0.3175  Median  Rank 

10 

0.290 

0.164 

0.082 

0.017 

C  =  0  Median  Rank 

10 

0.285 

0.220 

0.161 

0.080 

0.016 

C  =  0.375  Median  Rank 

10 

0.286 

0.221 

0.161 

0.081 

0.016 

C  =  0.3  (Standard)  Median  Rank 

10 

0.281 

0.212 

0.151 

0.076 

0.015 

C  =  1  Median  Rank 

20 

0.464 

0.392 

0.300 

0.176 

0.040 

Unmodified 

20 

0.468 

0.396 

0.303 

0.180 

0.040 

C  =  0.3175  Median  Rank 

20 

0.474 

0.401 

0.308 

0.183 

illCKl 

C  =  0  Median  Rank 

20 

0.3S5 

0.302 

0.178 

0.040 

C  =  0.375  Median  Rank 

20 

0.469 

0.396 

0.303 

wmm 

0.041 

C  =  0.3  (Standard)  Median  Rank 

20 

0.454 

0.381 

0.289 

0.170 

0.037 

0  =  1  Median  Rank 

30 

0.632 

0.553 

0.449 

0.299 

0.085 

Unmodified 

30 

0.637 

0.556 

0.451 

0.302 

0.087 

C  =  0.3175  Median  Rank 

30 

0.644 

0.562 

0.458 

0.307 

0.091 

C  =  0  Median  Rank 

30 

0.636 

0.554 

0.451 

0.301 

0.087 

C  =  0.375  Median  Rank 

30 

0.637 

0.556 

0.452 

0.302 

0.087 

C  =  0.3  (Standard)  Median  Rank 

30 

0.625 

tfism 

0.439 

0.289 

0.082 

C  =  1  Median  Rank 

40 

0.765 

0.698 

0.594 

0.431 

0.169 

Unmodified 

40 

0.767 

0.435 

0.172 

C  =  0.3175  Median  Rank 

40 

0.772 

0.707 

0.603 

0.445 

0.179 

C  =  0  Median  Rank 

40 

0.767 

0.700 

0.595 

0.433 

0.170 

C  =  0.375  Median  Rank 

40 

0.768 

0.701 

0.596 

0.436 

0.172 

C  =  0.3  (Standzurd)  Median  Rank 

40 

0.759 

0.690 

0.586 

0.424 

0.166 

C  =  1  Median  Rank 

50 

0.862 

0.808 

0.723 

0.564 

0.260 

Unmodified 

50 

0.864 

0.810 

0.726 

0.567 

0.263 

C  =  0.3175  Medizm  Rank 

0.867 

0.817 

0.729 

0.574 

0.271 

C  =  0  Median  Rank 

50 

0.863 

0.810 

0.726 

0.566 

0.262 

C  =  0.375  Median  Rank 

50 

0.864 

mm 

0.264 

C  =  0.3  (Standard)  Median  Rank 

50 

0.857 

0.800 

0.719 

0.552 

mi 

C  =  1  Median  Rank 

Table  4.14.  ANDERSON-DARLING  GOFT‘s  Power  Against  BETA(3,2)  data;  al 
pha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.219 

0.162 

0.108 

0.053 

0.012 

UnmodiOed 

10 

0.219 

0.163 

0.108 

0.053 

0.012 

C  =  0.3175  Median  Rank 

10 

0.220 

0.163 

0.110 

0.052 

0.013 

C  =  0  Median  Rank 

10 

0.219 

0.163 

0.108 

0.053 

0.012 

C  =  0.375  Median  Rank 

10 

0.219 

0.163 

0.109 

0.053 

0.012 

C  =  0.3  (Standard)  Median  Rank 

10 

0.216 

0.159 

0.106 

0.053 

\msi 

C  =  1  Median  Rank 

20 

0.256 

mm 

0.133 

0.072 

0.014 

Unmodified 

20 

0.258 

0.194 

0.134 

0.072 

0.014 

C  =  0.3175  Median  Rank 

20 

0.259 

0.197 

0.136 

0.073 

0.015 

C  =  0  Medizin  Rank 

20 

0.257 

0.194 

0.133 

0.072 

0.014 

C  =  0.375  Median  Rank 

20 

0.258 

0.194 

0.134 

0.072 

0.014 

C  =  0.3  (Standard)  Meaian  Rank 

20 

0.253 

0.186 

0.130 

0.070 

0.014 

C  =  1  Median  Rank 

30 

0.322 

0.249 

0.180 

0.101 

0.022 

Unmodified 

30 

0.327 

0.251 

0.182 

0.102 

0.023 

C  =  0.3175  Median  Rank 

lEm 

0.330 

0.253 

0.185 

0.103 

0.024 

C  =  0  Median  Rank 

0.326 

0.250 

0.182 

0.102 

0.023 

C  =  0.375  Median  Rank 

iKIil 

0.327 

0.251 

0.182 

0.102 

0.023 

C  =  0.3  (Standard)  Median  Rank 

30 

0.317 

0.243 

0.177 

0.098 

0.022 

C  =  i  Median  Rank 

40 

0.385 

0.318 

0.226 

0.132 

0.031 

Unmodified 

40 

0.386 

0.320 

0.228 

0.132 

C  =  0.3175  Median  Rank 

40 

0.388 

0.323 

0.232 

0.135 

0.034 

C  =  0  Median  Rank 

40 

0.386 

0.319 

0.227 

0.132 

0.031 

C  —  0.375  Median  Rank 

40 

0.386 

0.320 

0.228 

0.133 

0.032 

C  =  0.3  (Standeu-d)  Median  Rank 

\m\ 

0.381 

0.312 

0.223 

0.128 

rnrnsTiii 

0  =  1  Median  Rank 

n 

^1 

0.453 

0.373 

0.284 

0.163 

0.040 

Unmodified 

50 

0.456 

0.376 

0.287 

0.163 

0.040 

C  =  0.3175  Median  Rank 

50 

0.459 

0.379 

0.289 

mm\ 

0.042 

0  =  0  Median  Rank 

50 

0.455 

0.376 

0.286 

0.162 

mi 

C  =  0.375  Median  Rank 

50 

0.456 

0.376 

0.286 

0.163 

0.040 

C  =■  C.3  (Standard)  Median  Rank 

50 

0.447 

msMi 

0.282 

0.156 

0.036 

C  =  1  Median  Rank 
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Table  4.15.  ANDERSON-DARLING  GOFT's  Power  Against  EXPONENTIAL(l) 
data;  alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.668 

0.605 

0.525 

0.399 

0.215 

Unmodified 

10 

0.668 

0.604 

0.527 

0.399 

lEBiq 

C  =  0.3175  Median  Rank 

10 

0.667 

0.604 

0.526 

0.397 

C  =  0  Median  Rank 

10 

0.668 

0.604 

mm 

0.214 

C  =  0.375  Median  Rank 

10 

0.668 

0.604 

0.527 

0.399 

0.214 

C  =  0.3  (Standard)  Median  Rank 

10 

0.666 

0.603 

0.525 

0.399 

0.215 

C  =  1  Median  Rank 

20 

0.925 

0.902 

0.862 

0.783 

0.560 

Unmodified 

20 

0.925 

0.902 

0.863 

0.783 

0.561 

C  =  0.3175  Median  Rank 

20 

0.926 

0.903 

0.864 

0.783 

0.565 

C  =  0  Median  Rank 

20 

m 

0.903 

0.862 

0.784 

0.560 

C  =  0.375  Median  Rank 

20 

0.925 

0.902 

0.863 

0.783 

0.561 

C  =  0.3  (Standard)  Median  Rank 

20 

0.926 

0.901 

0.863 

0.782 

0.563 

C  =  1  Median  Rank 

30 

0.984 

0.978 

0.965 

0.935 

0.819 

Unmodified 

30 

0.985 

0.978 

0.965 

0.935 

0.819 

C  —  0.3175  Median  Rank 

30 

0.985 

0.978 

0.965 

0.934 

0.820 

C  =  0  Median  Rank 

30 

0.985 

0.978 

0.965 

0.935 

0.819 

C  =  0.375  Median  Rank 

30 

0.985 

0.978 

0.965 

0.935 

0.819 

C  =  0.3  (Standard)  Median  Rank 

30 

0.984 

0.977 

0.965 

0.936 

0.819 

C  =  1  Median  Rank 

40 

0.998 

0.996 

0.992 

0.983 

0.935 

Unmodified 

40 

0.998 

0.996 

0.992 

0.983 

0.935 

C  =  0.3175  Median  Rank 

40 

0.998 

0.996 

0.992 

0.983 

0.936 

C  =  0  Median  Rank 

40 

0.998 

0.996 

0.992 

0.983 

0.935 

C  =  0.375  Median  Rank 

40 

0.998 

0.996 

0.992 

0.983 

0.935 

C  =  0.3  (Standard)  Median  Rank 

40 

0.998 

0.995 

0.992 

0.983 

0.935 

C  =  1  Median  Rank 

50 

1.000 

1.000 

0.999 

0.995 

0.978 

Unmodified 

50 

1.000 

1.000 

0.999 

HI] 

0.979 

C  =  0.3175  Median  Rank 

50 

1.000 

1.000 

0.999 

im]m\ 

0.980 

C  =  0  Median  Rank 

50 

1.000 

1.000 

0.999 

0.995 

0.979 

C  =  0.375  Median  Rank 

50 

msM\ 

1.000 

0.999 

0.995 

0.979 

C  =  0.3  (Standard)  Median  Rank 

50 

1.000 

1.000 

0.999 

0.995 

0.978 

C  =  1  Median  Rank 
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Table  4.16.  ANDERSON-DARLING  GOFT's  Power  Against  DBL  EXPONEN 
TIAL(2,1)  data;  alpha  =  .20  ...  .01 


N 

.20  1 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

0.596 

0.485 

0.291 

Unmodified 

0.598 

0.485 

0.290 

C  =  0.3175  Median  Rank 

0.598 

0.482 

0.290 

C  =  0  Median  Rank 

0.598 

0.485 

0.291 

C  =  0.375  Median  Rank 

0.598 

0.485 

'  0.290 

C  =  0.3  (Standard)  Median  Rank 

0.597 

0.486 

0.290 

C  =  1  Median  Rank 

0.942 

0.911 

0.857 

0.687 

Unmodified 

0.942 

0.912 

0.858 

0.689 

C  =  0.3175  Median  Rank 

0.942 

0.911 

0.859 

0.692 

C  =  0  Median  Rank 

0.942 

0.911 

0.858 

0.688 

C  =  0.375  Mediam  Rank 

0.942 

0.912 

0.858 

0.689 

C  =  0.3  (Standard)  Median  Rank 

0.942 

0.913 

0.857 

0.688 

C  =  1  Median  Rank 

0.990 

0.984 

0.970 

0.990 

0.984 

0.970 

0.991 

0.984 

0.970 

0.990 

0.984 

0.970 

0.990 

0.984 

0.970 

0.990 

0.984 

0.969 

0.908  C  =  0.3175  Median  Rank 


0.908  C  =  0  Median  Rank 


0.908  C  =  0.375  Median  Rank 


1.000 

0.999 

0.997 

1.000 

0.999 

0.997 

1.000 

0.999 

0.997 

0.999 


0.999 


0.999 


0.997 


vmm 


0.993 


0.994 


0.993 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.997  0.993 


1.000 


1.000 


1.000 


1.000 


0.974 

Unmodified 

0.974 

C  =  0.3175  Median  Rani 

0.974 

C  =  0  Median  Rank 

0.974 

C  =  0.375  Median  Rank 

0.974 

C  =  0.3  (Standard)  Median  Rank 

0.974 

C  =  1  Median  Rank 

wiitini 


0.999 

0.994 

Unmodified 

0.999 

0.994 

C  =  0.3175  Median  Rank 

0.999 

0.994 

C  =  0  Median  Rank 

0.999 

0.994 

C  =  0.375  Median  Rank 

0.999 

0.994 

C  =  0.3  (Standard)  Medi2m  Rank 

0.999 

0.993 

C  =  1  Median  Rank  I 

Table  4.17.  ANDERSON-DARLING  GOFT's  Power  Against  WEIBULL(2,6)  data; 
alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.243 

0.190 

0.139 

0.075 

0.018 

Unmodified 

10 

0.243 

0.190 

0.140 

0.075 

0.018 

C  =  0.3175  Median  Rank 

10 

0.243 

0.189 

0.140 

0.074 

0.018 

C  =  0  Median  Rank 

10 

\mn 

0.189 

0.140 

0.075 

0.018 

C  =  0.375  Median  Rank 

10 

0.243 

0.190 

0.140 

0.075 

0.018 

C  =  0.3  (Standard)  Median  Rank 

10 

0.245 

0.190 

0.139 

0.076 

0.018 

C  =  )  Median  Rank 

20 

0.357 

0.288 

0.215 

0.131 

0.037 

Unmodified 

20 

0.357 

0.288 

0.214 

0.131 

0.037 

C  =  0.3175  Median  Rank 

20 

0.358 

0.288 

0.214 

fiireni 

0.037 

C  =  0  Median  Rank 

20 

0.357 

0.289 

0.214 

0.131 

0.037 

C  =  0.375  Median  Rank 

20 

0.357 

0.214 

0.131 

0.037 

C  =  0.3  (Standard)  Median  Rank 

20 

0.356 

0.284 

■ilHIil 

0.130 

0.037 

C  =  1  Median  Rank 

30 

0.454 

0.373 

0.284 

0.182 

0.058 

Unmodified 

30 

0.455 

0.373 

0.283 

0.182 

0.059 

C  =  0.3175  Median  Rank 

30 

0.457 

■ilCTa 

0.284 

0.182 

0.060 

C  =  0  Median  Rank 

30 

0.455 

0.372 

0.283 

0.182 

0.059 

C  =  0.375  Median  Rank 

30 

0.455 

0.373 

0.283 

0.182 

0.059 

C  =  0.3  (Standard)  Median  Rank 

30 

0.455 

0.368 

0.283 

0.180 

0.060 

C  =  1  Median  Rank 

40 

0.526 

0.454 

0.362 

0.241 

ItliCTlI 

Unmodified 

40 

0.526 

0.454 

0.363 

0.242 

0.091 

C  =  0.3175  Median  Rank 

40 

0.527 

0.453 

0.364 

0.244 

0.092 

C  =  0  Median  Rank 

40 

0.526 

0.454 

0.362 

0.241 

0.091 

C  =  0.375  Median  Rank 

40 

0.526 

0.454 

0.363 

0.242 

0.091 

C  —  0.3  (Standard)  Median  Rank 

40 

0.523 

0.453 

0.363 

0.240 

0.091 

C  =  1  Median  Rank 

50 

0.611 

0.531 

0.439 

0.305 

0.115 

Unmodified 

50 

0.612 

0.532 

0.440 

0.304 

0.116 

C  =  0.3175  Median  Rank 

1^1 

0.612 

0.532 

0.440 

0.305 

0.118 

C  =  0  Median  Rank 

50 

0.611 

0.531 

0.441 

0.304 

0.116 

C  =  0.375  Median  Rank 

50 

0.612 

0.532 

mi 

0.304 

0.116 

C  =  0.3  (Standard)  Median  Rank 

50 

0.609 

0.531 

0.439 

0.301 

0.111 

C  =  1  Median  Rank 
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Table  4.18.  ANDERSON-DARLING  GOFT's  Power  Against  LOG-NORMAL(0,1) 
data;  alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.779 

0.733 

0.677 

0.583 

0.406 

Unmodified 

10 

0.779 

0.732 

0.677 

im 

mm 

C  =  0.3175  Median  Rank 

10 

0.779 

0.730 

0.678 

0.582 

0.406 

0  =  0  Median  Rank 

10 

0.779 

0.732 

imn 

0.583 

0.406 

C  =  0.375  Median  Rank 

10 

0.779 

0.732 

1  0.677 

0.582 

B3 

C  =  0.3  (Standard)  Median  Rank 

10 

0.779 

0.732 

0.676 

0.583 

0.408 

C  =  1  Median  Rank 

20 

0.974 

0.963 

0.944 

0.905 

0.806 

Unmodified 

im 

0.975 

0.962 

0.944 

0.905 

0.806 

C  =  0.3175  Median  Rank 

Im 

0.975 

0.962 

0.944 

0.905 

0.808 

C  =  0  Median  Rank 

0.974 

0.963 

0.944 

0.905 

0.806 

C  =  0.375  Median  Rank 

1^ 

0.975 

iiiCTaa 

0.944 

0.905 

0.806 

C  =  0.3  (St2mdard)  Median  Rank 

20 

0.974 

0.962 

0.944 

0.905 

0.807 

C  =  1  Median  Rank 

Iglil 

0.998 

0.995 

0.992 

0.983 

0.950 

Unmodified 

30 

0.998 

0.995 

0.992 

0.983 

0.950 

C  =  0.3175  Median  Rank 

30 

0.998 

0.995 

0.992 

0.983 

0.950 

C  =  0  Median  Rank 

im 

0.998 

0.995 

0.992 

0.983 

0.950 

C  =  0.375  Median  Rank 

30 

0.998 

0.995 

0.992 

0.983 

0.950 

C  =  0.3  (Standard)  Median  Rank 

0.998 

0.995 

0.992 

0.983 

0.951 

C  =  1  Median  Rank 

□ 

1.000 

0.999 

0.998 

0.990 

Unmodified 

\m 

1.000 

1.000 

0.999 

0.998 

0.990 

C  =  0.3175  Median  Rank 

ESII 

1.000 

1.000 

0.999 

0.998 

0.991 

C  =  0  Median  Rank 

1.000 

1.000 

0.999 

0.998 

0.990 

C  =  0.375  Median  Rank 

mi 

1.000 

1.000 

0.999 

0.998 

0.990 

C  =  0.3  (Standard)  Median  Rank 

m 

1.000 

1.000 

0.999 

0.998 

0.990 

C  =  1  Mediem  Rank 

□ 

mi 

1.000 

1.000 

1.000 

1.000 

0.9S9 

Unmodified 

50 

1.000 

1.000 

1.000 

1.000 

0.999 

C  =  0.3175  Median  Rank 

50 

1.000 

1.000 

1.000 

ntililill 

0.999 

C  =  0  Median  Rank 

50 

1.000 

1.000 

1.000 

1.000 

0.999 

C  =  0.375  Median  Rank 

50 

1.000 

1.000 

1.000 

1.000 

0.999 

C  =  0.3  (Standard)  Median  Rank 

50 

1.000 

1.000 

1.000  1.000 

0.999 

C  =  1  Median  Rank 
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Table  4.19.  CRAMER- VON  MISES  GOFT's  Power  Against  UNIFORM(0,1)  data; 
alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.285 

0.226 

0.156 

0.079 

0.012 

Unmodified 

10 

0.306 

0.250 

0.175 

0.092 

0.015 

C  =  0.3175  Median  Rank 

10 

0.349 

0.285 

0.208 

0.115 

Mm 

C  =  0  Median  Rank 

10 

0.303 

0.242 

0.171 

0.088 

0.014 

C  =  0.375  Median  Rank 

10 

0.311 

0.252 

0.179 

0.093 

0.016 

C  =  0.3  (Standard)  Median  Rank 

m 

0.220 

0.164 

0.107 

IS^ 

0.008 

C  =  1  Median  Rank 

m 

0.426 

lip 

0.257 

0.146 

0.034 

Unmodified 

20 

0.459 

Illy 

0.285 

0.164 

0.039 

C  =  0.3175  Median  Rank 

20 

0.521 

0.438 

0.334 

0.204 

0.054 

C  =  0  Median  Rank 

20 

0.449 

0.373 

0.273 

0.158 

0.037 

C  =  0.375  Median  Rank 

20 

0.463 

0.384 

0.288 

0.168 

0.040 

C  =  0.3  (Standard)  Median  Rank 

20 

0.329 

0.257 

•>  ^74 

0.090 

0.018 

C  =  1  Median  Rank 

m 

0.552 

0.372 

0.214 

0.058 

Unmodified 

30 

0.585 

0.512 

0.406 

0.243 

0.066 

C  =  0.3175  Median  Rank 

30 

0.647 

0.566 

0.460 

0.295 

0.086 

C  =  0  Median  Rank 

30 

0.574 

0.499 

0.395 

0.234 

0.065 

C  =  0.375  Median  Rank 

30 

0.589 

0.517 

0.408 

0.245 

0.066 

C  =  0.3  (Standard)  Median  Rank 

30 

0.447 

0.366 

0.261 

0.147 

lililiTil 

C  =  1  Median  Rank 

40 

0.667 

0.584 

0.486 

0.330 

0.106 

Unmodified 

40 

0.696 

0.621 

0.523 

0.365 

0.118 

C  =  0.3175  Median  Rank 

40 

0.742 

0.678 

0.575 

0.413 

0.143 

C  =  0  Median  Rank 

40 

0.688 

0.611 

iiimi 

0.356 

0.115 

C  =  0.375  Median  Rjmk 

40 

0.699 

0.623 

0.526 

0.368 

0.120 

C  =  0.3  (Standard)  Median  Rank 

40 

0.569 

0.486 

0.382 

0.243 

0.063 

C  =  1  Median  Rank 

50 

0.768 

0.701 

0.603 

0.426 

0.176 

Unmodified 

50 

0.789 

0.728 

0.636 

0.460 

0.196 

C  =  0.3175  Median  Rank 

50 

0.826 

0.772 

0.686 

0.511 

0.237 

C  =  0  Median  Rank 

50 

0.783 

0.720 

0.626 

0.453 

0.191 

C  =  0.375  Median  Rank 

50 

0.790 

0.731 

0.638 

0.464 

0.198 

C  =  0.3  (Standard)  Median  Rank 

50 

0.686 

0.609 

0.503 

0.329 

0.120 

C  =  1  Median  Rank 
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Table  4.20.  CRAMER- VON  MISES  GOFT's  Power  Against  BETA(3,2)  data;  al 
pha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

0.176 

0.115 

0.056 

0.185 

0.123 

0.061 

0.199 

0.138 

0.068 

0.182 

0.120 

0.059 

0.186 

0.126 

0.061 

0.147 

0.097 

0.051 

Unmodified 


0.010  I  C  =  0.3175  Median  Rank 


0.011  C  =  0  Median  Rank 


C  =  0.375  Median  Rank 


0.260 

0.199 

0.130 

0.072 

0.015 

Unmodified 

0.271 

0.213 

0.139 

0.075 

0.017 

C  =  0.3175  Median  Rank 

0.301 

0.233 

0.160 

0.083 

0.021 

C  =  0  Median  Rank 

0.268 

0.208 

0.136 

0.074 

0.016 

C  =  0.375  Median  Rank 

0.273 

0.213 

0.139 

0.076 

0.017 

C  =  0.3  (Standard)  Median  Rank 

0.217 

0.159 

0.101 

0.057 

0.011 

C  =  1  Median  Rank 

0.239 

0.166 

0.089 

0.020 

Unmodified 

0.252 

0.178 

0.097 

C  s  0.3175  Median  Rank 

0.281 

0.202 

0.112 

0.024 

0  =  0  Median  Rank 

0.248 

0.175 

0.094 

0.022 

C  =  0.375  Median  Rank 

0.254 

0.179 

0.098 

0.021 

C  =  0.3  (Standard)  Median  Rank 

0.191 

0.131 

0.072 

0.015 

C  =  1  Median  Rank 

0.283 

0.208 

0.121 

0.303 

0.220 

0.130 

0.331 

0.250 

0.147 

0.298 

0.215 

0.128 

0.304 

0.222 

0.133 

0.237 

0.171 

0.093 

EMI 


0.030 


Unmodified 


.3175  Median  Rank 


C  =  0  Median  Rank 


C  =  0.375  Median  Rank 


0.428 


.461 


.423 


0.430 


.357 


0.336 

0.253 

0.137 

0.038 

Unmodified 

0.353 

0.266 

0.149 

0.042 

C  =  0.3175  Median  Rank 

0.382 

0.294 

0.167 

0.047 

C  =  0  Median  Rank 

0.347 

0.260 

0.145 

0.041 

C  =  0.375  Median  Rank 

0.355 

0.269 

0.150 

0.042 

C  =  0.3  (Standard)  Median  Rank 

0.288 

0.207 

0.108 

0.027 

0  =  1  Median  Rank 
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Table  4.21.  CRAMER- VON  MISES  GOFT‘s  Power  Against  EXPONENTIAL(l) 
data;  alpha  =  .20  ...  .01 


.15 

.10  i 

.05 

.01 

PLOTTING  POSITION 

10 

0.653 

0.593 

0.506 

0.383 

0.189 

Unmodified 

10 

0.652 

0.596 

0.509 

0.385 

0.182 

C  =  0.3175  Median  Rank 

10 

|i2j22l 

0.588 

0.507 

0.374 

0.174 

C  =  0  Median  Rank 

10 

0.655 

0.595 

0.507 

0.383 

0.183 

C  =  0.375  Median  Rank 

10 

0.596 

0.508 

0.384 

0.182 

C  =  0.3  (Standard)  Mediam  Rank 

10 

0.631 

0.576 

0.493 

0.378 

0.189 

C  =  1  Median  Rank 

iKtl 

0.903 

0.873 

0.818 

0.733 

0.519 

Unmodified 

20 

0.901 

0.870 

0.819 

0.726 

0.513 

C  =  0.3175  Median  Rank 

20 

0.897 

0.868 

0.815 

0.718 

0.503 

C  =  0  Median  Rank 

20 

wnmM 

0.871 

0.820 

0.728 

0.516 

C  =  0.375  Median  Rank 

20 

0.901 

0.870 

0.819 

0.726 

0.514 

C  =  0.3  (Standard)  Median  Rank 

20 

0.897 

0.865 

0.813 

0.729 

0.515 

C  =  1  Median  Rank 

0.971 


0.972 


0.971 


.971 


0.971 


0.962 

0.942 

0.895 

0.759 

Unmodified 

0.962 

0.942 

0.896 

0.754 

C  =  0.3175  Medi2in  Rank 

O.jH'  0 

0.943 

0.893 

0.745 

C  =  0  Median  Rank 

0.902 

0.942 

0.897 

0.759 

C  =  0.375  Medi2in  Rank 

0.961 

0.942 

0.895 

0.754 

C  =  0.3  (Standard)  Median  Rank 

0.959 

0.939 

0.891 

0.747 

C  =  1  Median  Rank 

0.993 


0.993 


0.993 


0.993 


.993 


0.989 

0.984 

0.964 

0.889 

Unmodified 

0.989 

0.965 

0.887 

C  =  0.3175  Medifin  Rank 

0.990 

0.985 

■ncisa 

0.883 

C  =  0  Median  Rank 

0.989 

0.984 

0.965 

0.888 

C  =  0.375  Median  Rank 

0.990 

itmi 

0.965 

0.887 

C  =  0.3  (Standard)  Median  Rank 

0.989 

0.982 

0.962 

0.886 

C  =  1  Median  Rank 

50 

0.999 

0.998 

0.996 

0.989 

0.959 

Unmodified 

50 

0.999 

0.998 

0.996 

0.989 

0.958 

C  =  0.3175  Median  Rank 

50 

0.999 

0.998 

0.989 

0.958 

C  =  0  Median  Rank 

50 

0.999 

0.998 

0.996 

0.989 

0.958 

C  =  0.375  Medicin  Rank 

50 

0.999 

0.998 

0.996 

0.989 

0.959 

C  =  0.3  (Standard)  Median  Rank 

50 

0.999 

0.998 

0.988 

0.956 

C  =  1  Median  Rank 
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fable  4.22.  CRAMER- VON  MISES  GOFT‘s  Power  Against  DBL  EXPONEN- 
TIAL(2,1)  data;  alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.708 

0.655 

0.577 

0.467 

0.254 

Unmodified 

10 

0.707 

0.656 

0.582 

0.469 

0.251 

C  =  0.3175  Median  Rank 

10 

0.702 

0.651 

0.574 

0.458 

0.240 

C  =  0  Median  Rank 

10 

0.708 

0.656 

0.580 

0.468 

113^ 

C  =  0.375  Median  Rank 

10 

0.708 

0.656 

0.582 

lESpS 

0.252 

C  =  0.3  (Standard)  Median  Rank 

10 

0.690 

0.637 

0.565 

0.256 

C  =  1  Median  Rank 

0.941 

0.920 

0.881 

0.819 

0.650 

Unmodified  1 

20 

0.941 

0.919 

0.881 

0.815 

0.643 

C  =  0.3175  Median  Rank 

20 

0.940 

0.916 

0.877 

0.806 

0.633 

C  =  0  Median  Rank 

20 

0.941 

0.920 

0.880 

0.817 

0.646 

C  =  0.375  Median  Rank 

20 

0.941 

0.919 

0.880 

0.815 

0.643 

C  =  0.3  (Standard)  Median  Rank 

20 

0.936 

0.916 

0.875 

0.813 

0.644 

0  =  1  Median  Rank : 

j 

30 

0.989 

0.981 

0.973 

0.947 

0.864 

Unmodified 

30 

0.989 

0.981 

0.973 

0.947 

0.861 

C  =  0.3175  Median  Rank 

30 

0.989 

0.981 

0.972 

0.946 

0.854 

C  =  0  Median  Rank  ' 

\m 

0.989 

0.981 

0.973 

0.947 

0.863 

C  =  0.375  Median  Rank 

30 

0.989 

0.981 

0.973 

0.861 

C  =  0.3  (Standard)  Median  Rank 

m 

0.987 

0.982 

lilRlil 

0.946 

0.859 

C  =  1  Medi£m  Rank  ! 

1  ■ 

0.998 

0.996 

0.994 

0.987 

0.951 

Unmodified  i 

\m 

0.998 

0.996 

0.994 

0.9o6 

0.950 

C  =  0.3175  Median  Rank 

40 

0.997 

0.996 

0.994 

0.986 

0.948 

C  =  0  Median  Rank 

40 

0.998 

0.996 

0.994 

0.987 

0.950 

C  =  0.375  Median  Rank 

40 

0.998 

0.996 

0.994 

0.986 

0.950 

C  =  0.3  (Standard)  Median  Rank 

40 

0.997 

0.996 

0.993 

0.986 

0.949 

C  =  1  Median  Rank 

50 

1.000 

1.000 

EEPI 

0.997 

0.989 

Unmodified 

50 

1.000 

1.000 

[iE^l 

0.997 

0.989 

C  =  0.3175  Median  Rank 

50 

um\ 

1.000 

mil 

0.997 

0.988 

C  =  0  Median  Rank 

50 

1.000 

1.000 

OB 

mail 

0.989 

C  =  0.375  Median  Rank 

50 

1.000 

1.000 

0.999 

0.997 

0.989 

C  =  0.3  (Standeird)  Median  Rank 

50 

1.000 

1.000 

1.000 

0.997 

0.987 

C  =  1  Median  Rank 
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Table  4.23.  CRAMER- VON  MISES  GOFT‘s  Power  Against  WEIBULL(2,6)  data; 
alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

10 

0.253 

0.203 

0.143 

0.072 

0.015 

Unmodified 

10 

0.257 

0.205 

0.145 

0.074 

0.015 

C  =  0.3175  Median  Rank 

10 

0.259 

0.208 

0.144 

0.077 

0.016 

C  =  0  Median  Rank 

10 

0.256 

0.204 

0.145 

0.073 

0.016 

C  =  0.375  Median  Rank 

10 

0.258 

0.204 

0.145 

0.074 

0.016 

C  =  0.3  (Standard)  Median  Rank 

10 

0.196 

0.135 

0.073 

0.015 

C  =  1  Median  Rank 

20 

0.339 

wma 

0.196 

0.115 

0.034 

Unmodified 

20 

0.338 

0.278 

0.199 

0.113 

0.033 

C  =  0.3175  Median  Rank 

20 

0.351 

0.278 

0.202 

0.113 

C  =  0  Median  Rank 

20 

0.338 

■iiwyil 

0.197 

0.114 

0.034 

C  =  0.375  Median  Rank 

20 

0.340 

0.200 

0.113 

0.033 

C  =  0.3  (Stcindard)  Median  Rank 

20 

0.331 

0.264 

0.191 

0.112 

0.034 

C  =  1  Median  Rank 

• 

m 

0.405 

0.334 

0.255 

0.154 

0.047 

Unmodified 

30 

0.410 

0.340 

0.256 

0.156 

0.047 

C  =  0.3175  Median  Rank 

30 

0.421 

0.344 

0.261 

0.158 

0.C46 

C  =  0  Median  Rank 

30 

0.409 

0.337 

0.256 

0.048 

C  =  0.375  Median  Rank 

30 

0.411 

0.340 

0.256 

0.155 

0.047 

C  =  0.3  (Standcird)  Mediam  Rank 

m 

0.395 

0.317 

0.237 

0.148 

0.043 

C  =  1  Mediem  Rank 

ictii 

0.468 

0.394 

0.316 

0.211 

0.066 

Unmodified 

40 

0.474 

0.402 

0.322 

0.212 

0.066 

C  =  0.3175  Median  Rank 

40 

0.484 

0.409 

0.324 

0.210 

C  =  0  Median  Rank 

40 

0.473 

0.401 

0.319 

0.212 

C  =  0.375  Median  Rank 

40 

0.475 

0.402 

inmMi 

0.213 

0.066 

C  =  0.3  (Standard)  Medijui  Rank 

40 

0.453 

0.383 

0.305 

0.202 

0.063 

C  =  1  Median  Rank 

50 

0.548 

0.475 

0.391 

0.254 

0.098 

Unmodified 

50 

0.554 

mil 

0.393 

0.255 

0.098 

C  =  0.3175  Median  Rank 

50 

0.563 

0.482 

0.399 

0.258 

0.099 

C  =  0  Median  Rank 

50 

0.554 

0.477 

0.391 

0.255 

0.099 

C  =  0.375  Median  Rank 

50 

0.555 

0.477 

0.394 

0.255 

0.099 

C  =  0.3  (Standard)  Medizm  Rank 

50 

0.535 

0.461 

0.373 

0.242 

0.096 

C  =  1  Median  Rank 

4-25 


Table  4.24.  CRAMER- VON  MISES  GOFT‘s  Power  Against  LOU-NORMAL(0,1) 
data;  alpha  =  .20  ...  .01 


N 

.20 

.15 

.10 

.05 

.01 

PLOTTING  POSITION 

0.768 


0.763 


.755 


0.766 


0.763 


0.963 


.963 


0.964 


0.963 


0.960 


0.727 

0.661 

0.567 

0.723 

0.662 

0.562 

0.711 

0.654 

0.546 

0.724 

0.662 

0.563 

0.722 

0.662 

0.561 

0.719 

0.659 

0.567 

0.369  I  Unmodified 


0.364  C  =  0.3175  Median  Rank 


C  =  0  Median  Rank 


0.367  I  C  =  0.375  Median  Rank 


0.363  I  C  =  0.3  (Standard)  Median  Rank 


0.948 


0.946 


0.948 


0.925 

0.883 

0.780 

Unmodified 

0.924 

0.881 

0.775 

C  =  0.3175  Medijui  Rank 

0.922 

0.877 

0.764 

C  =  0  Median  Rank 

0.924 

0.881 

0.777 

C  =  0.375  Median  Rank 

0.924 

0.881 

0.776 

C  =  0.3  (Standard)  Median  Rank 

0.922 

0.884 

0.779 

C  =  1  Median  Rank 

0.995 


0.995 


0.995 


0.995 


.995 


0.992 

0.987 

0.972 

Unmodified 

0.992 

0.987 

0.972 

0.926 

C  =  0.3175  Median  Rank 

0.991 

0.987 

0.971 

0.924 

C  =  0  Mediam  Rank 

0.99? 

0.987 

0.973 

0.929 

C  =  0.375  Mediim  Rank 

0.992 

0.972 

0.92f 

C  =  0.3  (Standard)  Medi2in  Rank 

0.991 

0.985 

0.973 

0.929 

C  =  1  Median  Rank 

0.999 


.999 


.999 


.999 


.999 


0.999 

0.998 

0.996 

0.982 

0.999 

0.998 

0.996 

0.981 

0.999 

0.999 

0.996 

0.981 

0.998 

0.996 

0  981 

0.999 

0.998 

0.996 

0.981  i 

0.999 

0.998 

0.996 

0.981 

C  =  0.3175  Median  Rank 


C  =  0  Median  Rank 


C  =  0.375  Median  Rank 


C  =  0.3  (Standard)  Median  Rank 


=  1  Median  Rank 


1.000 

1.000 

1.000 

0.999 

1.000 

WniiTil 

1.000 

0.999 

1.000 

1.000 

1.000 

0.999 

1.000 

1.000 

1.000 

G.999 

1.000 

1.000 

1.000 

0.999 

1.000 

1.000 

1.000 

0.999 

0.998 


0.998 


0.998 


Unmodified 


C  =  0.3175  Median  Rank 


C  =  0  Median  Rank 


C  =  0.375  Median  Rank 


C  =  0.3  (Standard)  Median  Rank 


C  =  1  Median  Rank 


V.  Conclusions 


This  chapter  is  a  presentation  of  the  author’s  conclusions  concerning  this  re¬ 
search  thesis. 

5.1  Review  of  Research  Objectives 

The  purpose  of  this  research  hcis  been  to  test  five  modifications  to  the  A-D  and 
the  C-VM  statistics  to  determine  if  an  increase  in  the  power  of  the  goodness-of-fit 
test  could  be  realized. 

5.2  Conclusions 

Two  genercJ  trends  are  apparent  from  the  power  tables  regardless  of  which  table 
is  being  examined.  The  first  is  that  as  the  sample  size  increases,  so  too  increases  the 
power.  This  makes  sense  because  as  more  information  is  gadned  about  a  population, 
it  becomes  more  apparent  that  it  is  not  from  the  normal  distribution.  The  second 
trend  is  that  as  a  increases  so  does  the  power  of  that  test.  The  reason  for  this 
is  not  quite  as  obvious,  but  to  shed  some  Jight,  consider  again  what  a  represents 
-  the  probability  of  wrongly  rejecting  the  null  hypothesis.  When  the  analyst  is 
less  concerned  about  erroneously  rejecting  a  true  hypothesis,  and  he  acts  on  this 
lessened  concern  by  selecting  a  larger  a,  then  it  follows  that  more  null  hypotheses 
will  be  rejected,  which  increases  the  overall  number  of  rejections,  and  subsequently, 
the  power! 

For  discussion  of  the  power  studies,  I  will  confine  my  analysis  primarily  to 
the  a  =  .05  levels  of  confidence  since  a  good  number  of  hypothesis  tests  use  this 
confidence  level.  One  side  note  concerning  the  validity  of  the  FORTRAN  subroutine 
POWERFUL  which  performs  the  power  studies  and  generates  the  tables;  when 
running  the  program  against  data  from  normad(0,l)  data,  as  expected,  the  power 
of  both  the  A-D  and  C-VM  tests  is  equal  to  a;  .20,  .15,  .10,  .05,  and  .01.  If  this 
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were  not  the  case,  then  something  would  be  wrong  because  a  is  the  probability  of 
rejecting  the  null  hypothesis  when  in  fcict  it  is  true. 

However,  when  testing  against  the  uniform(0,l)  distribution  ,  both  the  A-D 
and  the  C-VM  tests  are  fairly  powerful  with  the  A-D  being  the  stronger  of  the  two. 
The  strongest  individual  test  is  the  modified  A-D  test,  c=0,  n=50,  where  the  power 
is  .574.  At  all  a  and  n  levels,  the  modified  A-D  and  the  modified  C-VM  where  c=0 
are  the  most  powerful.  All  modifications  to  the  A-D  and  the  C-VM  statistics,  except 
for  c=l,  show  higher  power  than  the  tests  with  the  unmodified  statistics. 

It  is  a  different  story  for  the  test  powers  agsiinst  the  beta(3,2)  distribution. 
Neither  the  A-D  nor  the  C-VM  test  are  very  powerful.  The  strongest  individual  test 
is  the  modified  C-VM  test,  c=0,  n=50,  where  the  power  is  .167.  Again,  as  with  the 
uniform  distribution,  at  all  a  and  n  levels,  the  modified  A-D  and  the  modified  C-VM 
where  c=0  are  the  most  powerful.  All  modifications  to  the  A-D  statistic,  except  for 
c=l,  show  at  least  as  high  a  power,  if  not  higher  than  the  tests  with  the  unmodified 
A  D  statistic.  All  modifications  to  the  C-VM  statii.*ic,  except  for  c=l,  show  higher 
power  than  the  tests  with  the  unmodified  C-VM  statistic. 

Against  the  exponentiaJ(l)  distribution,  both  the  A-D  and  the  C-VM  tests  are 
very  powerful  with  the  A-D  being  ever-so-slightly  the  stronger  of  the  two.  There  is  a 
6-way  tie  for  the  strongest  individued  test.  It  is  amongst  all  of  the  A-D  tests  at  n=50 
where  the  power  is  .995.  Unlike  the  two  previous  distributions,  there  is  no  pattern  for 
any  one  statistic  dominating  the  power  studies  at  all  a  euid  n  levels,  parity  amongst 
the  statistics  seems  to  be  the  key  word  for  the  exponenti^J(l).  All  modifications  to 
the  A-D  and  the  C-VM  statistics  show  approximately  the  same  power  as  the  tests 
with  the  unmodified  statistics. 

Almost  identical  to  the  power  agfiinst  the  exponenti£d(l)  distribution,  for  the 
double  exponential(2,l,)  both  the  A-D  smd  the  C-VM  tests  are  very  powerful  with 
the  A-D  tests  being,  agciin,  slightly  stronger  than  the  C-VM  tests.  There  is  a  6-way 
tie  for  the  strongest  individual  test.  It  was  amongst  all  of  the  A-D  tests  at  n=50 


/ 
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where  the  power  is  .999.  Also  like  the  exponential(l),  there  is  no  app^Lrent  pattern  of 
any  one  statistic  dominating  the  power  studies  at  all  a  and  n  levels.  All  modifications 
to  the  A-D  and  the  C-V'M  statistics  show  approximately  the  same  power  as  the  tests 
with  the  unmodif  tatistics. 

The  po  .idy  against  the  Weibull(2,6)  distribution,  shows  both  the  A-D 
and  the  C-VM  tests  having  low  power.  For  the  small  sample  sizes,  there  is  a  rough 
similarity  between  the  powers  of  the  A-D  and  C-VM  tests,  but  as  the  sample  sizes 
approach  50,  the  A-D  test  is  clearly  the  stronger.  There  is  a  2-way  tie  for  the 
strongest  individual  test  between  the  unmodified  A-D  test  and  the  modified  A-D 
c=0  test,  at  n=50,  where  the  power  is  .305.  There  is  not  a  pattern  of  any  one 
statistic  dominating  the  power  studies  at  all  a  and  n  levels.  All  modifications  to  the 
A-D  and  the  C-VM  statistics,  except  for  c=l,  show  at  least  as  high  a  power  as  the 
tests  with  the  unmodified  statistics. 

Against  the  last  distribution  tested,  the  log-normal(0,l),  both  the  A-D  and 
the  C-VM  tests  are  very  powerful.  As  with  the  majority  of  the  other  distributions, 
the  A-D  test  is  slightly  more  powerful.  There  is  again  a  6-way  tie  for  the  strongest 
individual  test.  It  is  aunongst  adl  of  the  A-D  tests  at  n=50  where  the  power  is  .999. 
And  there  is  no  pattern  of  any  one  statistic  dominating  the  power  studies  at  adl  a  aind 
n  levels.  All  modifications  to  the  A-D  and  the  C-VM  statistics  show  approximately 
the  same  power  as  the  tests  with  the  unmodified  statistics. 

To  summaurize,  the  results  of  the  power  studies  against  the  exponential,  double 
exponentiad,  Weibull  and  log-normal  show  little  or  no  improvements  over  the  power 
of  the  unmodified  statistics  via  eight  of  the  modifications;  in  fact,  two  of  the  new 
statistics  show  a  marked  decrease  in  power  agaiinst  these  distributions.  In  the  case 
of  the  uniform  aind  beta  distributions,  the  power  increase  is  more  pronounced  with 
the  c=0  modification  to  the  A-D  and  C-VM  statistic  with  the  A-D  statistic  being 
the  stronger  (Overaill,  there  is  a  trend  of  the  A-D  tests  being  the  more  powerful 
versus  the  C-VM  tests).  Because  at  the  very  least,  no  power  is  lost  against  any  of 
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these  distributions,  but  actually  significantly  gciined  in  two  cases,  this  research  has 
shown  that  a  more  powerful  goodness-of-fit  test  for  normality  with  the  parameters 
estimated  from  the  sample  is  available  in  the  modified  A-D  statistic  defined  by 

ADmEDIAN,  C=0  =  -n-  2 
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Appendix  A.  FORTRAN  Programs 

«***4i4i««*«**«**«*«**«««*4i***«««4i«««*«**««*«**«««*****«4r4>«*«*«*« 

THIS  PROGRAM  DETERMINES  THE  CRITICAL  VALUES  WITH 
THE  MODIFIED  ANDERSON-DARLING  k  CRAMER-VON  MISES  STATISTICS 
VALUES  FOR  THE  NORMAL  DISTRIBUTION. 

THIS  PROGRAM  ALSO  GENiiRATES  THE  POWER  STUDIES. 

MANY  WRITE  STATEMENTS  HAVE  BEEN  COMMENTED  OUT  IN  THIS, 
CODE.  THE  PURPOSE  OF  THESE  WRITE  STATEMENTS  WAS  FOR  VERIF¬ 
ICATION  AND  TESTING  PURPOSES  AND  THEY  HAVVE  BEEN  LEFT  IN 
FOR  FUTURE  USERS  BENEFIT. 

MANY  THANKS  TO  KAHYA  AND  REAM  FOR  THE  BULK  OF  THIS  CODE. 
WITHOUT  THEIR  EFFORTS,  MY  THESIS  WOULD  NOT  HAVE  GONE  AS  FAR 
AS  IT  HAS. 

*******«**«««4««*«*****«*4i**4i«**«««  *«*««*«**««•****«*«**«***«*« 


C  •*  SOME  COMMENTS  ABOUT  THE  VARIABLES  USED  •• 

C 

C  S  -  •  OF  REPETITIONS,  10,000  THIS  THESIS 
C  N  -  SAMPLE  SIZE,  4,  5,  ...,  50 
C  ANDER( 10003)  -  AD  STATISTIC  VALUE  ARRAY 

C  SEED  ■  INITIAL  VALUE  NEEDED  TO  START  RANDOM  NUM  GENS 

THROUGH  IMSL  ROUTINE  RNSET  ( 

''XX(*)  -  OUTPUT  VECTOR  FOR  NORMAL  DEVIAT^ 

INTEGER  S,J,N,K,SS,SEED,III,I,RUNTYPE,PASSNO,IJ 
REAL  XX(50),ANDER(0;10003),DARL  ^ 

REAL  CRIT80 , CRIT85 , CRIT90 , CRIT95 , CRIT99 , AMORDF 
CHARACTER*32  COMBONAM 
EXTERNAL  RNNOR, RNSET, SVRGN,ANORDF 

C  THE  NUMBER  OF  REPETITIONS  IS 
S  -  10000 

PRINT  ♦ 

WRITE(*,1) 
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FORMAT ('THIS  PROGRAM  WILL  GEN  THE  AD  AND  CVM  GRIT  VALS  FOR  ONE', 

♦'  OF  THE  SEVEN  CHOICES  BELOW  (N  »  4  TO  50)  FOR  THE  N(0,1)  DISTR  ', 
♦'AND  THEN  PREPARE  A  POWER  STUDY  USING  DATA  FROM  THE  UNIF,  BETA  ', 
♦'  EXP,  DBL  EXP,  WEIB  AND  LNORM  DISTS  (N  ■  10,50,10).'//, 

♦'YOU  MUST  CHOOSE  FROM  THE  FOLLOWING  OPTS:',//, 

♦'1  -  USE  UNMODIFIED  PLOTTING  POSITION  WITH  AD  k  CVM  STATS',/, 

♦'2  -  USE  KAHYA'S  MESSED  UP  MEDIAN  RANK  PLOTTING  POSITION  ', 

♦'WITH  AD  k  CVM  STATS',/, 

♦'3  -  USE  C  -  0,3175  MEDIAN  RANK  PLOTTING  POSITION  * , 

♦'WITH  AD  *  CVM  STATS',/, 

♦'4  -  USE  C  -  0  MEDIAN  (KAHYA'S  MEAN)  RANK  PLOTTING  POSITION', 
♦'WITH  AD  k  CVM  STATS',/, 

♦'5  -  USE  C  -  0.375  MEDIAN  RANK  PLOTTING  POSITION  WITH  AD  ft  CVM  ' , 
♦'STATS',/, 

♦'6  -  USE  C  »  0.3  (STANDARD)  MEDIAN  RANK  PLOTTING  POSITION  WITH  ', 
♦'AD  ft  CVM  STATS',/, 

♦'7  -  USE  C  -  1  MEDIAN  RANK  PLOTTING  POSITION  WITH  AD  ft  CVM  STATS', 
♦//, 'ENTER  YOUR  SELECTION  AS  I .E.  3  OR  1 . '//) 

READ(*,2)RUNTYPE 
2  FORMAT(Il) 

IF(RUNTYPE.Eq.l) 

♦  COMBONAM-' UNMODIFIED  ' 

IF(RUNTYPE.Eq.2) 

♦  COMBONAM-' KAHYA'S  MESSED  UP  MEDIAN  RANK  ' 

IF(RUNTYPE.Eq.3) 

♦  COMBONAM-'C  -  0.3175  MEDIAN  RANK  ' 

IF(RUNTYPE.Eq.4) 

♦  COMBONAM-'C  •  0  MEDIAN  RANK  ' 

IF(RUNTYPE.Eq.5) 

♦  COMBONAM-'C  -  0.375  MEDIAN  RANK  ' 

IF(RUNTYPE.Eq.6) 

♦  COMBONAM-'C  -  0.3  (STANDARD)  MEDIAN  RANK' 

IF(RUNTYPE.Eq.7) 

♦  COMBONAM-'C  ■  1  MEDIAN  RANK  » 

IF (RUNTYPE . LT . 1 . OR . RUNTYPE . GT . 7)  STOP 
IF(RUNTYPE.Eq.l)THEN 

OPEN (UNIT- 1 , FILE- ' adcv 1 . dat ' , STATUS- ' UNKNOWN ' ) 

OPEN (UNIT-2 , FILE- ' c vmcvl . dat ' , STATUS- ' UNKNOWN ' ) 
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ENDIF 

IF (RUNTYPE . EQ . 2)THEN 

OPEN (UNIT- 1 , FILE- *  adcv2 . dat ' , STATUS- ' UNKNOWN ' ) 
OPEN  CUNIT-2 , FILE- ' cvmcv2 . dat ' , S  FATUS- ’ UNKNOWN ' ) 
ENDIF 

IF (RUNTYPE. EQ. 3) THEN 

OPEN (UNIT-1 , FILE- *  adcv3 . dat ' , STATUS- ' UNKNOWN ' ) 
OPEN (UNIT-2 . FILE- ’ cvmcvS . dat ' , STATUS- ' UNKNOWN ' ) 
ENDIF 

IF (RUNTYPE . EQ . 4) THEN 

OPEN (UNIT- 1 . FILE- *  adcv4 . dat ' , STATUS- ' UNKNOWN ’ ) 
OPEN (UNIT-2 , FILE- ’ cvmcv4 . dat ' , STATUS- ' UNKNOWN ' ) 
ENDIF 

IF (RUNTYPE. EQ. 5) THEN 

OPEN (UNIT- 1 , FILE- ’ adcvS . dat ' , STATUS- ' UNKNOWN ' ) 
OPEN (UNIT-2 , FILE- ' cvmcvS . dat ' , STATUS- ' UNKNOWN ' ) 
ENDIF 

IF (RUNTYPE . EQ . 6)THEN 

OPEN(UNIT-l , FILE- *  adcv6 . dat » , STATUS- ' UNKNOWN ' ) 
OPEN (UNIT-2 , FILE- » cvacvS . dat ' , STATUS- ' UNKNOWN ' ) 
ENDIF 

IF (RUNTYPE. EQ. 7) THEN 

OPEN (UNIT- 1 , FILE- ’ adcv7 . dat » , STATUS- ' UNKNOWN ' ) 
OPEN (UNIT-2 . FILE- » cvacv7 . dat ' , STATUS- ' UNKNOWN ' ) 
ENDIF 

DO  3000  PASSNO  -  1.2 

IF(PASSNO.Eq.l)THEN 
SEED  -  402958961 
ELSE 

SEED  -  510496597 
ENDIF 

CALL  RNSET(SEED) 

DO  2000  N  -  4,50 

•  INITIALIZE  EVERYTHING  FOR  THIS  RUN  OF  SAMPLE  SIZE  N.  • 


J-0 


A-3 


o  o  o  o  o  o  o 


K-0 

SS-0 

III-O 

I-O 

DO  5  I  -  l.N 
XX(I)-0.0 

5  CONTINUE 

DO  6  I  -  0,S+2 
ANDER(I}-0.0 

6  CONTINUE 
DARL-0.0 
CRIT80-0.0 
CRIT85-0.0 
CRIT90-0.0 
CRIT9S-0.0 
CRIT99-0.0 

0:tt******************************************** 

*  DO  THE  GRUNT  WORK  FOR  THIS  SAMPLE  OF  SIZE  N.* 

4i*«*4i**4i*4i4<4i4t****«i|i****«4i**«*#4i«**>ti**^*4i****4i*iti 

DO  11  J  -  l.S 
C 

C  CALL  RNNOR  TO  GENERATE  SET  OF  ORDER  STATISTICS 

C  FROM  THE  NORMAL  DISTRIBUTION 

C 

CALL  RNNOR (N. XX) 

CALL  STANDAR  TO  STANDARDIZE  THE  RANDOM  DEVIATES. 
THE  FOLLOWING  COMPUTATIONS  WILL 
USE  THESE  STANDARDIZED  VALUES. 

HRITE(*,*) 'NORMAL  DEVIATES  «  ' ,  (XX(IJ)  ,I>1 ,4) 
CALL  STANDAR (XX, N) 

WRITEC*,*) 'STANDARDIZED  DATA  -  ' , (XX(IJ) ,IJ-1 ,4) 

CALL  SVRGN  TO  SORT  THE  DATA  INTO  ASCENDING  ORDER 

CALL  SVRGN(N.XX.XX) 

C  WRITEC*,*) 'SORTED  DATA  ■  ',(XX(IJ),IJ-1,4) 
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DO  22  K-l.N 


C  USE  ANORDF  TO  DTRHN  THE  CDF  VALUES  OF  THE  NORMAL 

C  DISTRIBUTION  FOR  THE  STANDARDIZED  DATA. 

XX(K)  -  ANORDF(XX(K)) 

22  CONTINUE 

C  HRITE(*,*)*CDFED  DATA  (Uis)  ■  ',(XX(IJ),IJ-1,4) 


C  CALL  ANDERSON/CRAMER  TO  FIND  THE  AD  STAT  VALUES 

IF (PASSNO . EQ . 1) CALL  ANDERSON (N , XX , DARL , RUNTYPE) 

IF (PASSNO . EQ . 2) CALL  CRAMER (N , XX , DARL , RUNTYPE) 

ANDER(J)  -  DARL 

C  WRITE(*,*)'THIS  AD  STAT  IS  ',ANDER(J) 

11  CONTINUE 

ANDER(O)  -  0.0 
SS  •  S+1 

C  CALL  SVRGN  TO  SORT  THE  STATS  TO  DTRMN  THE  PERCENTILES 

C  OF  THE  ARRAY  IN  ORDER  TO  FIND  THE  CRITICAL  VALUES . 

CALL  SVRGN(SS,ANDER,ANDER) 

C  CALL  L;(TRA  to  XTRPLT  the  data  for  the  points  at  "0," 

C  AND  "N+1."  OUTPUT  VECTOR  IS  AN  ARRAY  OF  LENGTH  "N+2.'' 

CALL  EXTRA (S.ANDER) 

C  CALL  VALUES  TO  DETERMINE  THE  CRITICAL  VALUES 

CALL  VALUES (ANDER.CRITBO ,CRIT85 ,CRIT90 , CRIT95 , CRIT99 , SS) 
WRITE (PASSNO , 1 599 )  N , CRIT80 , CRIT85 , CRIT90 , CRIT95 , CRIT99 
1699  F0RMAT(1X,I2,5(2X,F6.4)) 

***************4r*****4i************ 

2000  CONTINUE 
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3000  CONTINUE 


DO  3010  1-1,2 
REVIND(I) 

3010  CONTINUE 

PRINT  • 

PRINT  The  critical  value  tables  for  the  menu 
'f'opt-on  you  selected  are  in  adcv(opt) .dat  ft  cvmcv(opt) .dat.  ' 
PR.N  ♦ 

CALL  POWERFUL (COMBONAM.RUNTYPE) 

PRINT  *,'The  Anderson-Darling  power  studies  for  the  uniform,  ', 
■•''beta,  exponential,  double  exponentizd  veibull,  and  lognormal  ', 
■•^'are  in  these  files:  adunif  (opt)  .dat,  adbeta(opt)  .dat,  ', 

■•' '  adexpn (opt )  .  dat ,  ' , 

•>’addblx(opt) .dat,  adweib(opt) .dat,  emd  adlogn(opt).dat. ' 

PRINT  * 

PRINT  *,'1110  Cramer- von  Mi ses  power  studies  for  the  uniform,  ', 
4'beta,  exponential,  double  exponential  weibull,  and  lognormal  *, 
<f*are  in  these  files:  cvmunif (opt) .dat,  cvmbota(opt) .dat,  ', 
'•''cvmexpn(opt)  .dat,  \ 

■f* cvmdblx(opt) .dat,  cvmveib(opt) .dat,  and  cvmlogn(opt) .dat. ' 

PRINT  ♦ 

STOP 

END 

* 

USING  THE  TECHNIQUE  EXPLAINED  IN  CHAPTER  3,  FIND  THE  * 
CRITICAL  VALUES  ♦ 

• 

SUBROUTINE  CAN(Y1,Y2,D1,D2,Y,RES) 

REAL  M,B,Y1,Y2,D1,D2,Y,RES 
IF((D2-Dl).Eq.0.0)D2  -  D2  ♦  1.00001 
M  -  (Y2-Y1)/(D2-D1) 
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B  -  Y1  “  M*D1 
RES  -  (Y-B)/M 
RETURN 
END 

*  * 

*  THE  FOLLOWING  SUB  CALCS  THE  AD  STAT  VALUES;  DEPENDING  ON  * 

*  THE  VALUE  OF  RUNTYPE,  THE  APPROPRIATE  MOD  IS  USED  (1:7).  ♦ 

*  * 
***m**^tti*:^mm^L^HL*^i^i**^^**iHL^i^^t^iti:^*it‘************************rHi*** 

SUBROUTINE  ANDERSON (K.R.ANDER, RUNTYPE) 

INTEGER  K, I, RUNTYPE 
REAL  R(50). TOTAL, XX,YY,ANDER,ZZ 
TOTAL  -  0.0 
c  print  * 

c  WRITE(*,100) 

c  print  ♦ 

c  100  FORMATC'N  I  I-l  Ui  In(Ui)  1-Un-i+l  ln(l-Un-i+l)  ' 
c  ♦»  ln(Ui)+ln(*, 

c  +'l-Un-i+l)  (I-l)*<ln(Ui)+ln(l-Un-i+l))  SUM'/, 


V.  -r  , 

C  +'  -  '/, 

c  ♦  '  N-1 

c  ♦'  N-1') 


DO  11  I  -  l.K 

IF(R(I).LE.O.O)  R(I)  «  0.0001 
ZZ  ■  1.0  -  R(K-I+1) 

IF(ZZ.LE.O.O)  ZZ  •  0.0001 
XX  -  LOGCRCD) 

YY  ■  LOG(ZZ) 

c  UNMODIFIED  PLOTTING  POSITION 

IF (RUNTYPE. EQ.l) 

♦  TOTAL  •  (2.0*REAL(I)-1.0)*(XX+YY)/(2.0*REAL(K))  +  TOTAL 

c  KAHYA'S  MESSED  UP  MEDIAN  RANK  PLOTTING  POSITION 

IF (RUNTYPE. EQ. 2) 

♦  TOTAL  -  (REAL(I)-0.3175)*(XX+YY)/(REAL(K)-0.365)  +  TOTAL 
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c  KAHYA'S  PROPER  MEDIAN  RANK  PLOTTING  POSITION 

IFCRUNTYPE.EQ.a) 

+  TOTAL  -  (REAL(I)-0.3175)*(XX+YY)/(REAL(K)+0.365)  +  TOTAL 

c  C  -  0  MEDIAN  (KAHYA'S  MEAN)  RANK  PLOTTING  POSITION 

IF(RUNTYPE.EQ.4) 

♦  TOTAL  •  (REAL(I)*(XX+YY))/(REAL(K)+1.0)  +  TOTAL 

c  GWINN'S  MEDIAN  RANK  PLOTTING  POSITION 

IFCRUNTYPE.EQ.S't 

+  TOTAL  -  (REAL(I)-0.375)*(XX+YY)/(REAL(K)+0.25)  +  TOTAL 

c  STANDARD  MEDIAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.6) 

+  TOTAL  »  (REAL(I)-0.3)*(XX+YY)/(REAL(K)+0.4)  +  TOTAL 

C  »  1  MEDIAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.EQ.7) 

TOTAL  -  (REAL(I)-1.0)*(XX+YY)/(REAL(K)-1.0)  +  TOTAL 
WRITE(*,150)K,I,(REAL(I)-1.0)/(REAL(K)-1.0),R(I),XX,ZZ,YY,XX+YY, 
(REAL(I)-1.0)*(XX+YY)/ (REAL (K ) - 1 . 0 ) , TOTAL 

F0RMAT(I1,2X,I1,2X,4(F7.4,1X),2X,F7.4,8X,F7.4,15X,F7.4,20X,F7.4) 

11  CONTINUE 
c  print  • 

c  print  *,'AD  -  -N  -2.  *  TOTAL’ 

c  PRINT  * 

ANDER  -  -REAL(K)  -  2.0  ♦  TOTAL 
RETURN 

END 

•*♦♦**♦♦***♦**♦*♦♦**♦*******♦******♦*♦♦♦***♦♦*•***♦*♦********♦*♦ 

*  * 

♦  THE  FOLLOWING  SUBR  CALCS  THE  CVM  STAT  VALUES;  DEPENDING  ON  • 

♦  THE  VALUE  OF  RUNTYPE,  THE  APPROPRIATE  MOD  IS  USED  (1:5).  ♦ 

♦  ♦ 


c 

+ 

c 

c  ♦ 
c  ISO 


SUBROUTINE  CRAMER(N,R,ANDER,RUNTYPE) 

INTEGER  I, N, RUNTYPE 

REAL  R(50), TOTAL, ANDER.MEDIAN 
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TOTAL  -  0.0 
DO  22  I  -  l.N 

c  UNMODIFIED  PLOTTING  POSITION 

IF(RUlfnVE.Eq.l) 

+  MEDIAN  -  (2.0*REAL(I)-1.0)/(2.0*REAL(N)) 

c  KAHYA'S  MESSED  UP  MEAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.2) 

♦  MEDIAN  «  (REAL(I)-0.3175)/(REAL(N)-0.365) 

c  KAHYA'S  PROPER  MEAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.3) 

♦  MEDIAN  -  (REAL(I)-0.3175)/(REAL(N)+0.365) 

c  C  -  0  MEDIAN  (KAHYA'S  MEAN)  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.4) 

♦  MEDIAN  -  REAL (I)/ (REAL (N) +1.0) 

c  GWINN'S  MEDIAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.5) 

+  MEDIAN  ■  (REAL(I)-0.375)/(REAL(N)+0.25) 

c  STANDARD  MEDIAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.6) 

+  MEDIAN  -  (REAL (I) -0.3)/ (REAL (N) +0.4) 

c  C  -  1  MEDIAN  RANK  PLOTTING  POSITION 

IF(RUNTYPE.Eq.7) 

♦  MEDIAN  -  (REAL(I)-1.0)/(REAL(N)-1.0) 

TOTAL  «  TOTAL  +  (R(I)  -  MEDIAN) **2 
22  CONTINUE 

ANDER  «  (1. 0/(12. 0*REAL(N)))  +  TOTAL 

RETURN 

END 

*  * 

*  THIS  SUBROUTINE  STANDARDIZES  ALL  THE  XX(I)  DATA  ♦ 

*  * 
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SUBROUTINE  STANDAR(X,N) 

INTEGER  I,N 

REAL  X(50),XSUM,XBAP  S.XOUT 
XSUM  -  0.0 
XOUT  -  0.0 
DO  100  I  -  l.N 
XSUM  -  XSUM  +  X(I) 

100  CONTINUE 

XBAR  «  XSUM/N 
DO  200  I  -  l.N 

XOUT  «  XOUT  +  (X(I)  -  XBAR)**2 
200  CONTINUE 

S  -  SQRT(XOUT/(N  -  1)) 

DO  300  I  -  l.N 
X(I)  -  (X(I)  -  XBAR)/S 
300  CONTINUE 
RETURN 
END 


*  * 

*  THIS  .SUBROUTINE  EXTRAPOLATES  THE  ANDER(I)  DATA  * 

•  TO  GENERATE  ANDER(O)  AND  ANDER(S+1)  FOR  COMPUTATION  ♦ 

♦  OF  THE  FIVE  CRITICAL  VALUES.  * 

♦  ♦ 


« 

SUBROUTINE  EXTRA (N.D) 

INTEGER  N.NO.Nl 

REAL  Y1.Y2.D(0: 10003) .D1.D2,ZZ 
Y1  -  0.5/|M 
Y2  -  1.5/N 
D1  ■  D(l)l 
D2  -  D(2) 1 

CALL  CAN(Y1.Y2.D1.D2.0.0.ZZ) 
IF(ZZ.GE.olo)  THEN 
D(0)  -  ZZ 
ELSE 

D(0)  -  0.0 
ENDIF  I 
Y1  -  (REAL(N)  -  1.5)/N 
Y2  -  (REAL (N)  -  0.5)/N 
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NO  -  N-1 
D1  -  D(NO) 

D2  -  D(N) 

CALL  CAN(Y1,Y2.D1,D2,1.0,ZZ) 

N1  -  N  +  1 

D(N1)  -  ZZ 

RETURN 

END 


*  * 

•  THE  FOLLOWING  SUB  DETERMINES  THE  XTILES  AND  FINDS  * 

•  THE  CRITICAL  VALUES  BY  EVOKING  THE  SUB  CAN  * 

•  * 


0:^*^^^^^**^i*^i:^i^i^:^**^**^Hi^^**^****m**^l^***********^^***************** 

SUBROUTINE  VALUES (D , CRIT80 , CRIT85 , CRIT90 , CRIT95 , CRIT99 , N) 
INTEGER  I,N,NN 

REAL  D(0 : 10003) ,Y(0 : 10003) ,C80 ,C90 ,C95 ,C99 .C85 , 

♦  Y79 , D79 , Y81 , D81 , DIF90 , Y89 , Y9 1 ,D89 , D91 , DIF95 , DIF80 , 

♦  Y94 , Y96 , D94 , D96 , DIF99 , Y98 , YlOO , D98 , DlOO , DIF85 , 

♦  Y84 , D84 , Y86 , D86 , CRIT85 , CRIT80 , CRIT90 , CRIT95 , CRIT99 
DO  100  I  «  1,N 

Y(I)  -  (REAL(I)  -  0.5)/REAL(N) 

100  CONTINUE 
Y(0)  -  0.0 
NN  -  N  +  1 
Y(NN)  -  1.0 
C80  -  1000.0 
C85  -  1000.0 
C90  -  1000.0 
C95  -  1000.0 
C99  -  1000.0 
DO  200  I  «  NN,0,-1 

IF  C'(I).LE.0.75)  GO  TO  300 
IF  (Y(I).GT.0.75.AND.Y(I).LE.0.80)  THEN 

C  GET  THE  DESIRED  XTILE  AT  807, 

DIF80  -  .80  -  Y(I) 

IF  (DIF80.LE.C80)  THEN 
C80  -  DIF80 


D79  -  D(I) 

Y81  -  Y(I+1) 

D81  -  D(I+1) 

ENDIF 

ELSEIF  (Y(I). GT. 0.80. AND. Y(I).LE. 0.85)  THEN 

GET  THE  DESIRED  XTILE  AT  857, 

DIF85  -  .85  -  Y(I) 

IF  (DIF85.LE.C85)  THEN 
C85  -  DIF85 
Y84  -  Y(I) 

D84  -  D(I) 

Y86  -  Y(I+1) 

D86  ■  D(I+1) 

ENDIF 

ELSEIF  (Y(I). GT. 0.85. AND. Y(I).LE. 0.90)  THEN 

GET  THE  DESIRED  XTILE  AT  90X 
DIF90  -  .90  -  Y(I) 

IF  (DIF90.LE.C90)  THEN 
C90  -  DIF90 
Y89  -  Y(I) 

D89  -  D(I) 

Y91  -  Y(I+1) 

D91  -  D(I+1) 

ENDIF 

ELSEIF  (Y(I).GT. 0.90. AND. Y(I).LE. 0.95)  THEN 

GET  THE  DESIRED  XTILE  AT  95% 

DIF95  -  .95  -  Y(I) 

IF  (DIF95.LE.C95)  THEN 
C95  -  DIF95 
Y94  -  Y(I) 

094  -  D(I) 

Y96  -  Y(I+1) 

D96  -  D(I+1) 

ENDIF 


ELSEIF  (Y(I). GT. 0.95. AND. Y(I).LE. 0.99)  THEN 


C  GET  THE  DESIRED  XTILE  AT  99y. 

DIF99  -  .99  -  Y(I) 

IF  (DIF99.LE.C99)  THEN 
C99  «  DIF99 
Y98  -  Y(I) 

D98  -  D(I) 

YlOO  -  Y(I+1) 

DlOO  -  D(I+1) 

ENDIF 

ENDIF 

200  CONTINUE 

300  IF  (DIF80.Eq.O.O)  THEN 
CRIT80  «  D79 
ELSE 

C  COMPUTE  THE  CRIT  VALUE  AT  SIGNIFICANCE  LEVEL  •  .20 

CALL  CAN(Y79,Y81,D79,D81, .80,CRIT80) 

ENDIF 

IF  (DIF85.EQ.0.0)  THEN 
CRIT85  ■  D84 
ELSE 

C  COMPUTE  THE  CRIT  VALUE  AT  SIGNIFICANCE  LEVEL  - 

CALL  CAN(Y84,Y86,D84,D86, .85,CRIT85) 

ENDIF 

IF  (DIFSO.EQ.O.O)  THEN 
CRIT90  -  D89 
ELSE 

C  COMPUTE  THE  CRIT  VALUE  AT  SIGNIFICANCE  LEVEL  -  . 10 

CALL  CAN(Y89,Y91,D89,D91, .90,CRIT90) 

ENDIF 

IF  (DIF95.Eq.0.0)  THEN 
CRIT95  ■  094 
ELSE 

C  COMPUTE  THE  CRIT  VALUE  AT  SIGNIFICANCE  LEVEL  ■  .05 

CALL  CAN(Y94, r96,D94,D96, .95,CRIT95) 

ENDIF 
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IF  (DIF99.Eq.0.0)  THEN 
CRIT99  -  D98 
ELSE 

C  COMPUTE  THE  CRIT  VALUE  AT  SIGNIFICANCE  LEVEL  -  .01 

CALL  CAN(Y98,Y100,D98,D100, .99,CRIT99) 

ENDIF 

RETURN 

END 


* 

THE  POWERFUL  SUB  EXECUTES  THE  POWER  STUDY  ♦ 
FOR  THE  ANDERSON-DARLING  AND  THE  CRAMER-VON  RISES  ♦ 
TESTS  WITH  THE  SIX  DISTRIBUTIONS  AT  DIFFERENT  ALPHA  ♦ 
LEVELS,  .20,  .15,  .10,  .05,  .01,  AND  AT  DIFFERENT  SAMPLE  * 
SIZES,  10,  20;  30,  40,  50  ♦ 


c  *****  SOME  COMMENTS  ABOUT  THE  VARIABLES  USED  •**♦* 

C  POWERC*)  -  THE  REJECTION  7, AGE  AT  THE  GIVEN  LEVEL 

SUBROUTINE  POWERFUL (COMBONAM,RUNTYPE) 

INTEGER  I , J , L , N , MAX , NUMDIST , RUNTYPE , IREPS IZE 
INTEGER  DISTRIB,STATTYPE, ALPHA, IREC 
INTEGER  SEED (6) 

REAL  S(50),CRITVAL(2,5,5),DUMMY(6) 

REAL  Y ,  P ,  ANDAR ,  AND ,  CRAMER  1 ,  CRAMER2 ,  POWER  ( (> ,  ?VS  ) ,  ANORDF 
CHARACTER*20  DISTNAM(6) 

CHARACTER*32  COMBONAM 
CHARACTER*34  STATNAM(2) 

EXTERNAL  RNUN,RNBET,RNEXP,RNEXT,RNWIB,SSCAL,RNLNL, 

*  SVRGN,ANORDF 

DATA  DISTNAM/' UNIFORM (0,1) » , 'BETA (3, 2) ' , ' EXPONENTIAL (1) ' 

+  ,'DBL  EXP0NENTIAL(2,1)','WEIBULL(2,6)','L0GN0RMAL(0,1)'/ 
DATA  STATNAM/' ANDERSON-DARLING  ', 

+  'CRAMER-VON  MISES  '/ 

DATA  SEED/77757985 , 172278979 , 536814143 , 751957587 , 4638907 , 

+  271411103/ 
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c 

c 


READ  1ST  RECORD  FROM  EA  FILE. 
THIS  CONTAINS  CRITVALS  FOR  N-4 


READ ( 1 .  •  )  (DUMK-f  (I )  ,  I-l  .6) 

READ(2,*)(DUHMY(I), 1-1,6) 

DO  10  STATTYPE  -  1,2 
DO  10  N  -  10,50,10 

C  SKIP  NEXT  FIVE  RECORDS 

DO  2  IREC  -  1,6 

READ (STATTYPE , *) (DUMMY (I ) , I - 1 , 6) 

2  CONTINUE 

READ(STATTYPE,*)DUMMY(1) , (CRITVAL (STATTYPE, ALPHA, N/ 10) , 
♦  ALPHA-1.5) 

C  SKIP  NEXT  FOUn  RECORDS 

IF(N.LT.50)THEN 
DO  7  IREC  •  1,4 

READ ( STATTYPE , • ) ( DUMMY ( I ) , I - 1 , 6 ) 

7  CONTINUE 

ENDIF 

10  CONTINUE 

IREPSIZE  -  5000 

MUMDIST  •  6 

DO  20  I  -  1,2 
CLOSE(I) 

20  CONTINUE 

MAX  -  50 

IF(RUNTYPE.Eq.l)THEN 

OPEN (UNIT-10 .FILE-' adunifl .dat ' .STATUS- 'unknown ' ) 

OPEN (UNIT- 1 1 , FILE- ' adbata 1 . dat ' , STATUS- ' unknown ' ) 

OPEN (UN IT- 1 2 , FI LE- ' adaxpn 1 . dat ' , STATUS - ' unknown ' ) 

OPEN (UNIT- 13 , FILE- ' addblx 1 . dat ' , STATUS- ' unknown ' ) 
0PEN(UNIT-14,FILE-'adw*ibl. dat '.STATUS- 'unknown') 
0PEN(UNIT-15,FILE-'adlognl. dat '.STATUS- 'unknown') 
0PEN(UNIT-20 .FILE- ' cvaunif 1 . dat ' , STATUS- ' unknown ' ) 

OPEN  (UNIT-2 1 ,  FILE- '  cviatat  al .  dat ' ,  STATUS- '  unknown ' ) 
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OPEN ( UN IT-2  2 , FI LE  ^ • c vmexpn 1 . dat ' , STATUS- ' unknown ' ) 
OPEN (UNIT-23 , FILE- ' cvndblxl . dat ' , STATUS- ' unknown ‘ ) 
OPEN (UNIT-24 , FILE- ' cvmwoibl . dat ' , STATUS- ' unknown ' ) 
OPEN (UNIT-25 .FILE- ' cvmlognl . dat ’ , STATUS- *  unknown ’ ) 
ENDIF 

IF(RUNTYPE.Eq.2)THEN 

OPEN (UNIT-10 .FILE- ' adunif 2 . dat ' .STATUS- 'unknown') 
OPEN (UN IT- 1 1 . FI LE- ' adb«ta2 . dat ' . STATUS- ' unknown ' ) 
0PEN(UNIT-12.FILE-'adexpn2.dat' .STATUS-'unknown') 
OPEN (UNIT-13 .FILE- ' addblx2 . dat ' . STATUS-'unknown ' ) 
0PEN(UNIT-14 .FILE- ' adweib2 . dat ' .STATUS- 'unknown ' ) 
OPEN (UNIT- 1 5 . FILE- ' adlogn2 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-20 . FILE- ' cvmunif 2 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-2 1 . FILE- ' cvmbet a2 . dat ' . STATUS- ' unknown ' ) 
OPEN  (UNIT-22  .FILE-'  cvjnexpn2 .  dat ' .  STATUS- '  unknown ' ) 
0PEN(UNIT-23 , FILE- ' cvindblx2 . dat ' . STATUS-'unknown ' ) 
OPEN  (UNIT-24  .FILE- '  cvinweib2 .  dat ' .  STATUS- '  unknown ' ) 
0PEN(UNIT-25 . FILE- ' cvBlogn2 . dat ' . STATUS- '«  nknown' ) 
ENDIF 

IF(RUNTYPE.EQ.3)THEN 

OPEN (UNIT- 1 0 . FILE- ' adunif 3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT- 1 1 . FILE- ' adb«ta3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT- 12 . FILE- ' adexpn3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-1 3 . FILE- ' addblx3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-14 .FILE- ' adweib3 . dat ' . STATUS-'unknown ' ) 
OPEN (UNIT- 1 5 . FILE- ' adlogn3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UN IT-20 , F I  LE- ' c VBuni f 3 . dat ' . STATUS- ' unknown  * ) 
OPEN  (UNIT-2 1 .  FILE- '  cvinb«ta3 .  dat ' .  STATUS- '  unknown ' ) 
OPEN (UNIT-22 .FILE- ' cvraexpn3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-23 .FILE- ' cvmdblx3 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-24 . FILE- ' cv«weib3 . dat ' . STATUS- ' unknown ' ) 
OPEN  rUNIT-qS . FILE- ' cvmlogn3 . dat ' . STATUS-'unknown ' ) 
ENDIF 

IF(RUNTYPE.Eq.4)THEN 

OPEN (UNIT- 10 .FILE- ' adunif 4 . dat ' , STATUS- 'unknown ' ) 
0PEN(UNIT-1 1 .FILE- ’ adb«ta4 . dat ' , STATUS- ' unknown ' ) 
OPEN  (UNIT- 12 .  FILE- '  ad#‘xpn4 .  dat ' .  STATUS- '  unknown ' ) 
OPEN (UNIT-13 .FILE- ' addl \x4 . dat ' . STATUS- 'unknown ' ) 
0PEN(UNIT-14.FILE-'adw«ib4. dat '.STATUS- 'unknown') 
OPEN (UNIT- 1 5 .FILE- ' adlogn4 . dat ' . STATUS- ' unknown ' ) 
0PEN(UNIT-20 .FILE- ' cvmunif 4 . dat ' , STATUS- ' unknown' ) 
OPEN (UNIT-2 1 , FILE- ' cvmb«ta4 . dat ' , STATUS- ' unknown ' ) 


OPEN (UNlT-22 .FILE- * cvmexpn4 . dat* , STATUS  ' unknown ' ) 
OPEN (UNlT-23 . FILE- *  cvndblx4 . dat ' . STATUS- ’ unknown ' ) 
OPEN (UNIT-24 .FILE- *  cvaweiM . dat ' , STATOS- ' unknown ' ) 
OPEN (UNIT-25 .FILE- *  cvmlogn4 . dat ' . STATUS- ' unknown ' ) 

ENDIF 

IF(RUNTYPE.Eq.5)THEN 

OPEN (UN IT- 1 0 . FILE- *  adun if  5 . dat ' , STATUS- ' unknown ' ) 
OPEN (UNIT- 1 1 , FILE- *  adbetaS . dat ' , STATUS- ’ unknown ' ) 
0PEN(UNIT-12 . FILE- ’ adexpnS . dat  * . STATUS- 'unknown ’ ) 
0PEN(UNIT-13 .FILE- ' addblxS . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-14 . FILE- ' adwoi b5 . dat ' . STATUS- ' unknown ' ) 
0PEN(UNIT-15,FILE-'adlogn5. dat '.STATUS- 'unknown') 
OPEN (UNIT-20 .FILE- ' cvmunif 5 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-2 1 .FILE- ' cvnbetaS . dat ' . STATUS- ' unknown ' ) 
L.  t2l (UNIT-22 , FILE- ' cvnexpnS . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-23 , FILE- ' cvndblxB . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-24 , FILE- ' cvaweibS . dat ' . STATUS- ’ unknown ' ) 
0PEN(UNIT-25 .FILE- ' cvalognB . dat ' , STATUS- ' unknown ' ) 

ENDIF 

IF(RUNTYPE.Eq.6)l’HEN 

OPEN (UNIT-10 .FILE- 'adunif 6 . dat ' .STATUS- 'unknown ' ) 
OPEN (UNIT-1 1 . FILE- ' adbetaS . dat ' . STATUS- ' unknown ' ) 
0PEN(UNIT-12 .FILE-' adexpnS . dat ' .STATUS- 'unknown ' ) 
OPEN (UNIT-13 . FILE- ' addblxS ,  'at ' . STATUS- ' unknown ' ) 
0PEN(UNIT-14.FILE-'adweibS.  *  .STATUS- 'unknown') 
0PEN(UNIT-15.FILE-'adlognS.da; » .STAPJS- 'unknown') 
OPEN (UNIT-20 . FILE- ' cvaunif S . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-2 1 .FILE- ' cvabetaS . dat ' , STATUS- ' unknown ' ) 
OPEN (UNIT-22 . FILE- ' cvnexpnS . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-23 . FILE- ' cvndblxS . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-24 . FILE- ' cvnweibS . dat ' . STATUS- ' unknown ' ) 
OPEN  (UNIT-25 .  FILE- '  cvnlognS .  dat ' ,  ST  *.TUS- '  unknown ' ) 

ENDIF 

lP(RUNTYPE.Eq.7)THEN 

0PEN(UNIT-10.FILE-'adunif7. dat '.STATU '-'unknown') 
OPEN (UNIT- 1 1 , FILE- ' adbeta7 . dat ' . STATUl  • ' unknown ' ) 
OPEN (UNIT- 12 . FILE- ' adexpn7 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT- 13 . FILE- ' addblx7 . dat ’ , STATUS- ' unknown ' ) 
0PEN(UNIT-14.FILE-'adweib7. dat '.STATUS- 'unknown') 
OPEN (UNIT- 1 6 . FILE- ' adlogn7 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-20 . FILE- ' cvmunif 7 . dat ' . STATUS- ' unknown ' ) 
OPEN (UNIT-2 1 , FILE- ' cvmbeta7 . dat ' , STATUS- ' unknown ' ) 
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OPEN (UNIT-22 .FILE- *  cvnexpn? . dat ' , STATUS- 'unknown ' ) 
OPEN (UNIT-23 , FILE- ' cvmdblxT . dat ' , STATUS- ' unknown ' ) 
OPEN (UNIT-24 .FILE- ' cvmweibT . dat ' .STATUS- 'unknown' ) 
OPEN (UNIT-25 .FILE- ' cvalogn? . dat ' . STATUS- 'unknown ' ) 
ENOIF 

DO  1000  N  -  10. MAX. 10 

•  INITIALIZE  EVERYTHING  FOR  THIS  RUN  OF  SAMPLE  SIZE  N.  * 

****«*««**4>*******4i*i»**«**«***«*4i«**4>***«**««*****4>****«* 


DO  60  DISTRIB  -  l.NUMDIST 
DO  60  STATTYPE  -  1.2 
DO  50  ALPHA  -1.5 

POWER(DISTRIB. STATTYPE. ALPHA)  -  0.0 

50  CONTINUE 

Y-0.0 

P-0.0 

ANDAR-0.0 

AND-0.0 

CRAMERl-0.0 

CRAMER2-0.0 

DO  51  I  -  l.N 
S(I)-0.0 
61  CONTINUE 

DO  600  DISTRIB  -  l.NUMDIST 
CALL  RNSET(SEED(DISTRIB)) 

DO  600  J  -  I.IFJIPSIZE 

C  GENERATE  THE  RANDOM  DEVIATES! II!!! 

GO  T0(60.65.70,76,80.85)DISTRIB 


o  o 


^  WRITE(*.55) 

55  FORMAT('E3mOR  AT  GENERATION  COMPUTED  GO  TO. ') 

STOP 

C  UNIFORM, BETA, EXPON.DOUBL  EXPON,WEIBULL,  ft  LOGNORMAL 

C  DETAILS  ON  THESE  RANDOM  i  GENERATORS  CAN  BE  HAD  FROM 

/  C  THE  IMSL  STAT/LIBRARY  MANUAL  WITH  THE  RESPECTIVE  PAGES 

C  P963,  P993,  P999,  PlOOl,  P1025,  PlOll. 

60  CALL  RNUN(N,S) 

GO  TO  90 

65  CALL  RNBETCN, 3. 0,2.0, S) 

GO  TO  90 

70  CALL  RNEXP(N,S) 

GO  TO  90 

75  CALL  RNEXT(N,2.0,1.0,0.5,S) 

GO  TO  90 

80  CALL  RNVIB(N,2.0.S) 

CALL  SSCAL(N,6.0,S,1) 

GO  TO  90 

85  CALL  RNLNL(N,0.0,1.0,S) 

C  STANDARDIZE  AND  SORT  THE  STANDARDIZED  DATA! II ! I! ! 

95  F0RMAT(A15,2X,I2,/,10(F5.2,1X)) 

90  CALL  STAN0AR(S,N} 

CALL  SVRGH(N,S,S) 

DO  400  L  -  1,N 
Y  -  SCL) 

DETERMINE  THE  CDF  VALUES 
FOR  THE  STANDARDIZED  DATA. 

S(L)  -  ANORDF(S(L)) 

400  CONTINUE 

C  FIND  THE  STAT  VALUES  FOR  THE  FOLLOWING  TESTS. 

CALL  ANDERSON (N,S,ANDAR,RUNTyPE) 
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CALL  CRAMER (N.S.CRAMERl.RUNTYPE) 


C  ENTER  THE  CRITICAL  VALUES  FOR  THE  SAMPLE  OF  SIZE  K  ’.'HM 

C  THE  CRITICAL  VALUE  TABLES. 

DO  450  ALPHA  -  1,5 

IF  (ANDAR  . GT. CRITVALd, ALPHA, N/ 10)) 

♦  P0WER(DISTRIB,1, ALPHA)  -  P0WER(DISTRIB,1, ALPHA)  +1.0 

IF  (CRAMER1.GT.CRITVAL(2, ALPHA, N/10)) 

+  P0HER(DISTRIB,2,ALPHA)  -  P0WER(DISTRIB,2,ALPHA)  +  1.0 

450  CONTINUE 


500  CONTINUE 

*******^*«**<***4i*««****<*4>4'4>***4>*4i*4i**4>**«*«****4t**«««*******4i*«*i»* 

•  .  I  .  * 

•  DETERMINE  THE  REJECTION  XAGE,  THE  POWER  * 

•  !  * 
0**0*0t*********************************0************************* 


DO  550  STATTYPE  -  1,2  j 

DO  550  ALPHA  -  1,5  | 

POWER(DISTRIB , STATTYPE , ALPHA) -POWERCDISTRIB , 

♦  STATTYPE, ALPHA)/REAL(IREPSIZE) 

650  CONTINUE  ! 


600  CONTINUE 

•  ♦ 

•  PRINT  THE  HEADERS  AT  THE  TOP  OF  ALL  12  FILES.  • 

*  * 


IF(N.EQ.10)THEN 
DO  606  fiTATTYPE  -  1,2 
DO  60L;  DISTRIB  -  1,NUMDIST 
WRITE(STATTYPE*10+DISTRIB-1 ,603)DISTNAM(DISTRIB) , 
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•f 

603 

♦ 

♦ 

♦ 

♦ 

605 


STATNAM(STATTYPE)  .COMBOHAl! 

FORMAT (4X, 'ACTUAL  POPULATION:  ',1(20.11, 
2X, 'STATISTIC:  ',A34,//. 

2X, 'PLOT  POSITION  IS  ',A32,///, 
,7X, 'POWERS  AT  ALPHA  LEVELS',//, 
'N',3X,'.20  .15  .10  .05 

CONTINUE 

ENDIF 


.01'/) 


^m^^t^^nm********************************************************** 


* 

* 

* 


PRINT  THE  POWERS 


^1*0**^************************************************************ 


DO  620  STATTYPE  -  1,2 
DO  620  DISTRIB  -  1,NUMDIST 
C  PRINT  THE  REJECTION  XAGE,  THE  POWER 

WRITE (STATTYPE* 10+DISTRIB- 1 , 6 10) N , (POWER (DISTRIB , STATTYPE , 
*  ALPHA), ALPHA-1, 5) 

610  F0RMAT(I2,2X,5(F5.3,3X)) 

620  CONTINUE 

1000  CONTINUE 

DO  1100  1-1,2 
DO  1100  J  -  1,6 
CL0SE(I*10+J-1) 

1100  CONTINUE 

RETURN 

END 

*  * 

*  THE  END  OF  THE  ROAD! M ! I !  « 

♦  .  * 
(I***«******«*«**4i4'*«4t*«*«*««******««*««****«**4i^«4i*«4i  4,41 
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